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ABSTRACT 


This  investigation  attempts  to  experimental ly  determine  the 
horizontal  shear  stress  distribution,  from  the  longitudinal  strain 
distribution,  across  the  depth  of  a  glulam  beam.  The  maximum  mid-depth 
experimental  shear  stress  values  are  then  compared  to  theoretical  mid¬ 
depth  shear  stresses  in  the  vicinity  of  a  mid-depth  delamination. 

One  5"  x  30  7/8"  x  40'  -  0"  Douglas  Fir  glulam  beam  was 
extensively  instrumented  and  tested,  within  the  elastic  range,  under 
twenty-eight  different  test  conditions.  Varying  conditions  included; 
delamination  length,  span/depth  ratio,  and  position  of  the  delamination 
with  respect  to  moment  and  shear. 

Beam  tests  indicated  that  material  properties,  load  proximity 
and  loading  conditions  effected  the  flexural  stress  distributions  and 
hence  also  effected  the  shear  stress  distributions. 

Theoretical  shear  stresses  were  found  to  be  tolerably  accurate 
at  any  section  a  distance  greater  than  one  quarter  the  beam  depth  from 
the  end  of  a  mid-depth  delamination.  On  sections  closer  to  the  mid-depth 
delamination  inadeauate  delamination  length,  a  short  gauge  length  and 
load  proximity  did  not  permit  a  comparison  between  the  experimental  and 
theoretical  mid-depth  shear  stresses. 

A  series  of  tests  are  recommended  to  check  the  magnitude  of  the 
shear  stresses  within  one  quarter  the  beam  depth  from  the  end  of  a  mid¬ 
depth  delamination. 
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CHAPTER  I 


INTRODUCTION 


A.  Introductory  Remarks 

In  recent  years  one  topic  of  interest  in  glued  laminated 
timber  beam  research  has  been  the  effect  of  the  stress  concentration 
in  the  region  of  a  delamination.  In  1960  Tamberg  and  Huggins  (13) 
concluded  from  a  small  number  of  tests  that  a  delamination  could 
seriously  reduce  the  strength  of  a  glulam  beam.  Huggins,  Aplin  and 
Palmer  (5)  concluded  from  futher  tests  that  shear  stresses  in  the 
vicinity  of  a  mid-depth  delamination  could  be  computed  with  "tolerable 
accuracy".  Collins  (2)  investigated  the  computation  of  the  shear 
stresses  by  experimentally  defining  a  strain  profile. 

The  major  drawback  to  all  of  these  previous  tests,  however, 
is  that  the  test  specimens  were  so  shallow  that  even  minor  defects  or 
minor  experimental  errors  played  a  significant  role  in  the  apparent 
behaviour  of  the  beam.  It  was  for  this  reason  that  the  present  study 
of  the  behaviour  of  a  full-size  glulam  beam  was  initiated. 

B.  Object 

The  object  of  the  present  program  was  to  determine  the  long¬ 
itudinal  shear  stresses  in  the  vicinity  of  a  delamination. 

The  longitudinal  shear  stress  may  be  derived  from  flexural 
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strains.  Longitudinal  strain  profiles  at  two  sections  in  a  beam  span 
may  be  converted  into  flexural  stress  distributions  by  applying  the 
stress  -  strain  relationship  for  the  material.  From  the  resulting 
stress  profiles  shear  stresses  between  the  two  sections  may  then  be 
calculated  at  various  locations  throughout  the  beam  depth. 

In  the  test  program  the  following  variables  were  examined: 

1.  delamination  length 

2.  span/depth  ratio 

3.  position  of  delamination  with  respect  to  moment  and  shear 
C.  Scope 

A  single  beam  specimen  was  tested  under  twenty-eight  differ¬ 
ent  test  conditions.  The  beam  was  composed  of  19  -  1  5/8  inch 
laminations  and  had  a  finished  size  of  5"  x  30  7/8"  x  40'-  0".  Tests 
were  conducted  for  a  single  load  at  the  mid-span  and  for  two  loads 
symmetrical ly  placed  about  the  mid-span  of  the  beam.  Tests  were 
carried  out  for  no  delamination  and  for  delamination  lengths  of  0.5, 
1.0,  and  1.5  times  the  beam  depth.  Span/depth  ratios  of  8,  10,  and 
12  were  investigated. 
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CHAPTER  II 


REVIEW  OF  PREVIOUS  RESEARCH 

A.  Experimental 

Prehaps  the  first  research  on  the  subject  of  shear  strength 
of  checked  or  split  timber  beams  was  conducted  by  Newlin,  Heck  and 
March  (9)  in  1951.  In  an  experimental  investigation  of  200  art¬ 
ificially  checked  timber  beams  it  was  noted  that  the  usual  method  of 
calculating  longitudinal  shearing  stresses  in  the  neutral  plane  did 
not  accurately  predict  the  ultimate  shearing  stresses.  It  was 
concluded  that  the  upper  and  lower  halves  of  the  beam  acted  to  some 
extent  independently  and  hence  relieved  some  shearing  stress  on  the 
neutral  plane. 

In  1959  Longworth  and  Stieda  (8)  investigated  the  effect  of 
a  delamination  in  a  region  of  a  scarf  joint.  It  was  found  that  the 
delamination  tended  to  lower  both  the  flexural  strength  and  the  stiff¬ 
ness  of  the  beam. 

In  a  test  program  involving  70  beams  in  1960,  Tamberg  and 
Huggins  (13)  concluded  that  delamination  resulted  in  a  significant 
reduction  in  the  load  carrying  capacity  of  the  beam. 

Strong  (12)  carried  out  an  extensive  experimental  program 
on  239  beams  in  1964,  with  position  and  length  of  the  delamination 
as  major  variables.  It  was  concluded  that  the  length  and  position 
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of  the  delamination  are  significant  factors  in  the  load  carrying 
capacity  of  a  beam;  end  delamination  is  more  critical  than  interior 
delamination  and  a  mid-depth  delamination  is  more  critical  than 
delamination  at  1/4  or  3/4  depth  when  the  delamination  is  at  the  mid¬ 
span  of  the  beam.  It  was  shown  that  variations  in  the  specific  gravity 
of  individual  laminations  effect  the  stress  distribution. 

Huggins ,Apl in  and  Palmer  (5)  concluded  from  tests  on  127 
beams,  that  an  average  shear  stress,  calculated  on  the  undelaminated 
surface,  did  not  explain  the  relative  incidence  of  horizontal  shear 
failures. 

Collins  (2),  in  1966,  attempted  to  calculate  horizontal 
shear  stress  from  a  strain  profile  obtained  by  instrumenting  3"  x  6" 
x  6'  -  2"  beams  at  uniform  intervals  over  the  beam  length.  Shear 
stresses  were  calculated  on  the  basis  of  constant  modulus  of  elasti¬ 
city.  It  was  concluded  however  that  the  beams  tested  were  not 
sufficiently  deep  to  accurately  define  the  strain  profile. 

B.  Theoretical 

A  method  of  evaluating  shear  stresses  in  delaminated  beams 
was  developed  by  Schwaighofer ,  Lee  and  Microys  (11)  and  is  based 
upon  work  reported  by  Iyengar  (6).  The  procedure  is  divided  into 
two  individual  calculations  (FIGURE  2-1). 

The  first  step  in  the  solution  is  to  determine  the  flexural 
stress  distribution.  This  is  done  by  writing  two  simultaneous 
equations  in  terms  of  the  moment  and  the  compressive  force  in  the 
portion  of  the  beam  above  the  delamination.  The  first  equation  is 
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FIGURE  2.1  SCHEMATIC  DIAGRAM  OF  THEORETICAL  SHEAR  STRESS  CALCULATION 
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written  using  the  static  moment  equilibrium  while  the  second  equation  is 
obtained  by  considering  the  strain  distribution  in  the  top  and  bottom 
portions  of  the  beam.  After  the  unknowns  have  been  determined,  the  flexural 
stress  distribution  can  be  solved  by  using  the  standard  flexural  formula. 

The  second  step  in  the  solution  is  the  determination  of  the  shear 
stresses.  This  portion  of  the  problem  is  based  on  a  method  developed  by 
Guyon  (4)  who  approximated  shear  stresses  along  a  post-tensioned  prestressed 
concrete  beam  by  multiplying  the  magnitude  of  the  prestressing  force  by  an 
influence  coefficient.  This  influence  coefficient  is  dependent  on  the 
position  of  the  prestressing  force  with  respect  to  the  location  where  the 
shear  stresses  are  to  be  determined. 

Similarly  shear  stresses  in  a  delaminated  beam  may  be  calculated 
by  employing  influence  coefficients.  To  date  influence  coefficients 
developed  by  Guyon  (4)  and  by  Lee  (7),  (see  FIGURE  2-2),  have  been  used 
to  predict  the  mid-depth  shear  stresses  in  glulam  beams.  To  determine  the 
shear  stress  on  an  elemental  length  the  area  of  the  flexural  stress  distri¬ 
bution  diagram  is  multiplied  by  the  area  of  the  appropriate  influence  line 
diagram. 

Huggins,  Aplin  and  Palmer  (5)  compared  the  influence  coefficient 
proposed  by  Guyon  (4)  and  by  Lee  (7)  and  concluded  that  Lee's  influence 
coefficients  predicted  the  shear  stresses  with  "tolerable  accuracy". 


ffKrtJod  bnB  qo*  9rt:t  n"  rv  ?n  dHiarb  rrbnJa  sril  pnh:>breno3  \  ‘  i  snr&Jdo 

* 


rte  bn.  '  sil  <  ivscn 


- fi ‘ u  6*  V  bn.  bfts  e  srtt  pnuu  \o  be*',  Fo;  9d  nep  no  Jud  ip?  b  aae'  Ja 


b9229n^29*tq  bonorzneJ-Jzoq  6  onofs  2S2rnie  n  9rla  bsismi xonc;qs  ortw  (£)  no^bt) 


sdj  9P9i  m  norPepoT  9riP  oP  }D9q29T  riPFw  sonoP  pnfzaynJapnq  9rtP  "to  norPraoq 

. 

■ 


INFLUENCE  COEFFICIENT  "C 


7 


SHEAR  STRESS 


CP 

2B 


FIGURE  2.2  INFLUENCE  LINES  FOR  SHEAR  STRESSES 


CHAPTER  III 


TEST  BEAM  SPECIMEN 


A.  Introduction 

In  order  to  obtain  a  well  defined  flexural  strain  profile  it 
was  necessary  to  have  a  deep  beam  section.  For  obvious  economical 
reasons  the  number  of  specimens  to  be  tested  had  to  be  kept  to  a 
minimum.  As  a  result  it  was  decided  to  use  one  large  beam  and  to  test 
this  beam  under  a  variety  of  conditions. 

B.  Material  Selection 

The  laminating  stock  used  was  1  5/8"  x  5  1/2"  Douglas  Fir 
classified  as  Canadian  Standard  Association  24f  grade.  The  stock  was 
sorted  with  respect  to  moisture  content,  approximate  specific  gravity, 
grain  orientation,  and  occurrence  of  knots  and  checks. 

The  selected  material  was  check  free,  exceptionally  free  of 
knots  and  generally  flat  grained.  Material  was  selected  on  the  basis 
of  a  moisture  content  of  10%  ±  2%.  The  approximate  specific  gravities 
of  individual  pieces,  determined  by  weighing,  ranged  from  a  maximum  of 
0.62  to  a  minimum  of  0.40  with  the  majority  of  the  stock  in  the  0.46  to 
0.54  range. 

C.  Fabrication  of  the  Beam 

Standard  fabrication  procedures  were  employed.  The  beam  was 
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laid  up  with  the  highest  specific  gravity  in  the  outer  lamination  decreasing 
to  the  lowest  specific  gravity  at  mid-depth.  Scarf  joints  occurred  at 
random  throughout  the  beam.  No  joints  occurred  in  the  area  of  the  delami¬ 
nation.  The  glulam  beam  was  finished  to  a  dimension  of  5"  x  30  7/8"  x  40' 

-  0"  (FIGURE  3-1). 

Casein  glue  at  72°F  was  applied  at  the  rate  of  62  lbs./  1000  sq„ 
ft..  During  gluing  operations  the  room  temperature  was  68°F  and  the 
relative  humidity  was  58%.  Open  assembly  time  for  the  beam  was  14  minutes. 
An  assembly  pressure  of  100  psi  was  applied  for  a  period  of  24  hours  (PLATE 
3-1). 
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FIGURE  3.1  BEAM  SPECIMEN 


CHAPTER  IV 


DESCRIPTION  OF  BEAM  TESTS 


A.  Introduction 

The  major  part  of  the  program  consisted  of  testing  the  beam 
with  a  delamination  located  at  mid-span  under  varying  conditions  of 
delamination  length,  span/depth  ratio,  bending  moment  and  shear. 

The  program  was  planned  to  include  additional  tests  with  the  delamin¬ 
ation  located  at  the  quarter  point  of  the  span  under  the  same  varying 
conditions  as  the  mid-span  delamination  series.  To  provide  for  these 
additional  tests,  using  the  same  glulam  beam,  the  following  test 
arrangement  was  devised. 

A  beam  1.25  times  the  length  of  the  longest  span  required 
for  the  mid-span  delamination  series  was  fabricated.  By  locating 
the  delamination  2/5  of  the  length  from  one  end  of  the  beam,  tests 
could  be  conducted  with  a  delamination  located  at  the  quarter  point 
of  the  span  or  at  the  mid-span.  The  different  test  conditions  are 
shown  schematically  in  FIGURE  4-1. 

To  balance  the  weight  of  the  overhanging  end  of  the  beam 
a  counterweight  was  placed  on  the  beam.  The  position  of  the 
counterweight  did  not  vary  for  different  span/depth  ratios  for  a 
given  series. 
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C.L.  LOAD 

C.L.  DELAMINATION 


LOAD  TYPE  1 


2  PT.  LOAD 
C.L.  DELAMINATION 


LOAD  TYPE  2 


1/4  PT.  LOAD 

1/4  PT.  DELAMINATION 


LOAD  TYPE  3 


C.L.  LOAD 

1/4  PT.  DELAMINATION 


LOAD  TYPE  4 


3/4  PT.  LOAD 

1/4  PT.  DELAMINATION 


LOAD  TYPE  5 


FIGURE  4.1  LOADING  CONDITIONS 
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B.  Loading  System 

Loads  were  applied  by  means  of  a  standard  testing  frame  and 
the  "Amsler  Hydraulic  Jacking  System"  available  in  the  Graduate 
Structural  Laboratory  (FIGURE  4-2  and  PLATE  4-1A).  Load  was  trans¬ 
ferred  to  the  test  beam  from  the  jack  through  a  6"  x  1"  x  1'  -  0" 
distribution  plate.  Reactions  were  transmitted  to  the  beam  through 
6"  x  1"  x  1'  -  0"  bearing  plates  which  were  free  to  rotate  in  a 
vertical  plane  (  PLATE  4-2). 

Lateral  support  for  the  beam  was  provided  at  two  locations. 
A  roller  system  was  clamped  to  both  sides  and  to  the  top  and  bottom 
of  the  beam  (  PLATE  4-3).  The  rollers  were  supported  by  vertical 
channels  bolted  to  the  testing  frame.  Movement  of  the  beam  was 
thereby  restricted  to  the  vertical  plane.  Friction  was  kept  to  a 
minimum  by  lubricating  the  rollers  and  the  flanges  of  the  channels. 

C.  Instrumentation 

i.  Deflection  Measurement 

Vertical  deflections  were  obtained  from  readings  on  scales 
suspended  from  the  beam  at  five  uniformly  spaced  positions. 
Deflections  were  determined  to  an  accuracy  of  ±  0.005"  by  means  of 
a  "Sokkisha  Novus  Automatic  Level". 

ii.  Longitudinal  Strain  Measurement 

An  8-inch  Demec  gauge  was  modified  to  read  strains  on  a 
gauge  length  of  3.90  inches.  The  modification  was  made  by  clamping 
a  conversion  bar  securely  to  the  Demec  gauge  (  PLATE  4-1  B). 
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FIGURE  4.2  TYPICAL  BEAM  TESTING  ARRANGEMENT 


A.  TYPICAL  TEST  SETUP 


B.  CONVERTED  DEMEC  GAUGE 


PLATE  4.1 
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B.  ROLLER  SUPPORT 


PLATE  4.2 


A. 


CLOSE  UP  OF  LATERAL  SUPPORT 


B.  OVERALL  VIEW  OF  LATERAL  SUPPORTS 


PLATE  4.3 
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The  calibration  factor  for  the  converted  Demec  strain  gauge  was  obtained 
from  a  tension  test  of  a  1/2  inch  diameter  steel  rod.  The  rod  was 
instrumented  with  two  electrical  resistance  SR-4  Type  A-3  strain  gauges 
and  two  sets  of  Demec  gauge  points.  Deformations  on  the  3.90  inch  gauge 
length  were  determined  to  an  accuracy  of  ±  0.00001  in/in. 

Demec  points  were  bonded  to  the  test  beam  with  sealing  wax  (FIGURE 
4-3).  The  number  of  gauge  lines  employed  in  each  test  was  dependent  on 
the  length  of  the  delamination.  For  tests  involving  no  delamination  and  a 
delamination  length  of  1.5  times  the  depth  of  the  beam  strains  on  gauge 
lines  1  through  30  were  read.  For  the  delamination  length  of  0.5  times  the 
beam  depth,  strains  were  read  for  gauge  lines  1  through  16  and  26  through  30. 
Strains  were  read  for  gauge  lines  1  through  20  and  26  through  30  for  a  delami¬ 
nation  length  equal  to  1.0  times  the  beam  depth. 

FIGURE  4-4  shows  the  position  of  the  mid-depth  delamination  within 
the  instrumented  section  of  the  beam. 

D.  Test  Conditions 

TABLE  4-1  indentifies  the  tests  conducted. 

The  maximum  moment  applied  to  the  beam  was  the  same  for  all  values 
of  span/depth  ratio  and  for  all  loading  conditions.  The  Canadian  Standards 
Association  Standard  086-1959  allows  an  extreme  fiber  bending  stress  of  2400 
psi  on  C.S.A.  24f  grade  material.  Load  increments,  based  on  C.S.A. 
specifications,  were  chosen  to  give  extreme  fiber  stresses  of  800  psi,  1600 
psi,  2400  psi  and  3200  psi. 
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FIGURE  4.3  BEAM  INSTRUMENTATION 
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A.  DELAMINATION  =  0.5  DEPTH 


GAUGE  LINES 
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B.  DELAMINATION  =  1.0  DEPTH 
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C.  DELAMINATION  =1.5  DEPTH 
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FIGURE  4.4  POSITION  OF  MID-DEPTH  DELAMINATION 


L/D  DELAMINATION  LENGTH  DELAMINATION  LOAD 

RATIO  NONE  0.5D  1.0D  1.5D  LOCATION  CONDITION 
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E.  Testing  Procedure 

Initial  Demec  strain  gauge  readings  were  recorded  the  evening 
previous  to  the  day  the  test  was  conducted. 

Load  increments  were  applied  at  a  rate  of  1000  pounds  per  minute. 
Deflection  readings  and  Demec  strain  gauge  readings  were  recorded.  A 
waiting  period  of  approximately  20  minutes  was  reguired  before  strain  gauge 
readings  were  commenced  for  the  last  load  increments.  Duration  of  the  beam 
tests,  involving  four  load  increments,  varied  from  5  to  7  hours. 

The  tests  were  conducted  in  such  a  sequence  that  shear  conditions 
became  more  severe  as  the  program  progressed.  For  all  delamination  lengths 
the  two  point  load  test  preceded  the  single  point  load  test  for  each  span/ 
depth  ratio.  Span/depth  ratios  were  tested  in  order  from  the  highest  to 
the  lowest. 

On  completion  of  the  first  group  of  tests,  involving  no  delami¬ 
nation,  the  beam  was  removed  from  the  test  frame  and  a  delamination  was 
cut  by  means  of  a  circular  power  saw  and  squared  off  with  a  hand  saw. 

After  the  delamination  was  cut  Demec  points  were  replaced  because  of 
excessive  wear.  Further  increases  in  delamination  length,  when  required, 
were  made  with  a  hand  saw.  Demec  points  were  again  replaced  before  the 
start  of  the  last  series  of  tests  in  which  the  delamination  length  was 
1.5  times  the  beam  depth. 

F.  Recording  of  Strain  Data 

Demec  strain  gauge  readings  were  recorded  directly  on  previously 
indentified  computer  mark  sense  cards. 
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At  each  load  increment  a  minimum  of  440  to  a  maximum  of  600 


strain  readings  were  recorded  depending  on  the  length  of  delamination. 


CHAPTER  V 


MATERIAL  PROPERTIES 

A.  Introduction 

On  completion  of  the  test  program  a  5'  -  0"  length  was  cut  from 
the  instrumented  section  of  the  beam.  This  section  was  then  cut  into 
the  19  laminations.  Each  5'  -  0"  length  lamination  was  then  cut  into 
lengths  of  3'  -  0"  and  2'  -  0". 

B.  Modulus  of  Elasticity  in  Compression 

Compression  test  specimens  were  fabricated  from  the  2'  -  0" 
lengths.  Each  length  was  ripped  in  half  and  each  half  was  planed  to 
a  uniform  thickness  of  1-inch.  The  two  corresponding  pieces  were  glued 
together  with  casein  glue  under  a  pressure  of  100  psi  applied  for 
twenty-four  hours,  and  then  trimmed  to  final  dimension  of  2"  x  2"  x  1'  - 
2"  (FIGURE  5-1  A). 

The  compression  test  specimens  were  loaded  to  failure  in  a 
"Baldwin-Lima-Hamil ton  Testing  Machine"  of  300,000  pound  capacity 
(PLATE  5-1  A).  Longitudinal  strains  over  an  8-inch  gauge  length  were 
obtained  by  means  of  a  "H.F.  Moore  Mechanical  Extensometer" .  The 
determination  of  the  modulus  of  elasticity  is  outlined  in  Appendix  A. 


25 


!  !  O',  H  [  i-  '  i.qfr  <J  '1 


r  -  ■  ■  -  •'  ;-f 


.(a  i-a  3ainn)  us 


26 


irtl  = 

> 

CN 

u3  -  » » 1 02 

Csl 

CO 


TENSION  SPECIMEN 


A.  COMPRESSION  SPECIMEN  TEST 


B.  TENSION  SPECIMEN  TEST 


PLATE  5.1 


C.  Modulus  of  Elasticity  in  Tension 


Tension  specimens  were  fabricated  from  the  3'  -  0"  lengths. 
Each  length  was  ripped  in  half  and  one  half  was  reduced  to  an  uniform 
thickness  of  1/2  inch.  The  1/2  inch  piece  was  shaped  and  trimmed  to 
the  finished  size  (FIGURE  5-1  B).  2  1/2"  x  1/2"  x  6"  pieces  cut  from 

stock  were  then  glued  to  the  ends  of  the  shaped  specimen  with  casein 
glue  and  held  under  pressure  for  twenty-four  hours. 

The  tension  test  specimens  were  loaded  to  failure  in  a 
"Bal dwin-Tate-Emery  Universal  Testing  Machine"  of  200,000  pound  capa¬ 
city  (PLATE  5-1  B).  Longitudinal  strains  were  measured  by  means  of 
Demec  gauge  with  an  8-inch  gauge  length  (see  Appendix  A). 

D.  Moisture  Content  and  Specific  Gravity 

Moisture  content  and  specific  gravity  was  determined  for 
each  lamination  on  the  basis  of  tests  performed  on  1  1/4"  x  1  1/4" 
x  1/2"  samples.  The  methods  employed  are  outlined  in  Appendix  B. 


CHAPTER  VI 


PRESENTATION  OF  RESULTS 

A.  Material  Properties 

The  test  results  of  the  laminating  stock  are  presented  in  TABLE 
6-1,  FIGURES  6-1,  6-2  and  6-3  and  PLATE  6-1. 

B.  Test  Results 

The  beam  test  results  are  presented  as  follows: 

1.  Variations  in  delamination  length  (PLATE  6-2) 

2.  Maximum  deflections  based  on  load  equal  to  1.33  times  the  working 
load  (TABLES  6-2,  6-3  and  6-4) 

3.  The  effect  of  load  proximity  on  strain  distribution  (FIGURE  6-4 
4„  Typical  strain  distributions  not  influenced  by  load  proximity 

(TABLE  6-5  and  FIGURE  6-5  through  FIGURE  6-9) 

5.  The  effect  of  the  modulus  of  elasticity  on  equilibrium  conditions 
(TABLE  6-6) 

6.  An  example  of  shear  stress  distribution  (FIGURES  6-10  and  6-11) 

7.  Mid-depth  shear  stresses  (FIGURE  6-12  through  FIGURE  6-20) 

8.  Effect  of  combined  shear  and  bending  moment  on  equilibrium 
conditions  (TABLE  6-7) 

9.  Horizontal  shear  failure  (FIGURES  6-21  and  6-22  and  PLATES  6-3 
and  6-4) 
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A  computer  program  was  written  to  determine  the  shear 
at  the  mid-depth  of  the  beam.  A  description  of  this  program  i 
Appendix  C. 


stresses 
found  in 


A  complete  record  of  data  and  the  detailed  results  have  been 
placed  on  file  with  the  Department  of  Civil  Engineer!" no. 


INITIAL  FINAL  SPECIFIC  ULTIMATE  STRESS  MODULUS  OF  ELASTICITY 

MOISTURE  MOISTURE  GRAVITY  TENSION  COMPRESSION  TENSION  COMPRESSION 

CONTENT  CONTENT  (KSD  (KSI)  (KSI)  (KSI) 
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TABLE  6.1  MATERIAL  PROPERTIES 
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DEPTH  LENGTH  TEST  BEND  AVERAGE  MODIFIED 
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TABLE  6.6  MODULUS  OF  ELASTICITY  EFFECT  ON  EQUILIBRIUM 
CONDITIONS  IN  A  PURE  MOMENT  REGION 


SPAN  DELAMINATION 

DEPTH  LENGTH  TWO  POINT  LOAD  MID-SPAN  LOAD 
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TABLE  6.7  EQUILIBRIUM  CONDITION  COMPARISON  BETWEEN  A  REGION  OF  PURE 
MOMENT  AND  A  REGION  OF  COMBINED  MOMENT  AND  SHEAR 
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FIGURE  6.1  VARIATION  OF  MOISTURE  CONTENT  AND  ULTIMATE  STRENGTH  ON  CROSS-SECTION  OF  BEAM 
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FIGURE  6.2  VARIATION  OF  SPECIFIC  GRAVITY  AND  MODULUS  OF  ELASTICITY  ON  CROSS-SECTION  OF  BEAM 
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FIGURE  6.3  GRAPH  OF  MODULUS  OF  ELASTICITY  VS  SPECIFIC 
GRAVITY  OF  INDIVIDUAL  LAMINATIONS 
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FIGURE  6.4  EFFECT  OF  LOAD  PROXIMINITY  ON  STRAIN  PROFILES 
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FIGURE  6.5  TYPICAL  STRAIN  DISTRIBUTION  FOR  ALL  INCREMENTS  OF  LOAD 
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A.  PURE  MOMENT 


SPAN/DEPTH  RATIO  =12  O  O  DEMEC  STRAINS 

GAUGE' LINE  11  -  FITTED  STRAINS 


FIGURE  6.6  WORKING  LOAD  STRAIN  DISTRIBUTION 

NO  DELAMINATION 
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FIGURE  6.7  WORKING  LOAD  STRAIN  DISTRIBUTION 

DELAMINATION  LENGTH  =  0.5D 
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FIGURE  6.8  WORKING  LOAD  STRAIN  DISTRIBUTION 

DELAMINATION  LENGTH  =  1.0D 


B.  COMBINED  MOMENT  AND  SHEAR 
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FIGURE  6.9  WORKING  LOAD  STRAIN  DISTRIBUTION 

DELAMINATION  LENGTH  =  1.5D 
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A.  STRAIN  DISTRIBUTION 


GAUGE  LINE  20 
GAUGE  LINE  16 


100 

B.  FITTED  STRESS  DISTRIBUTION 


SPAN/DEPTH  RATIO  =  10  WORKING  LOAD 

NO  DE LAM I NAT I ON 


FIGURE  6.10  AN  EXAMPLE  OF  STRAIN  AND  STRESS  DISTRIBUTION 

IN  A  REGION  OF  COMBINED  MOMENT  AND  SHEAR 


A.  NON-CORRECTED  SHEAR  STRESS  PROFILE 
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THEORETICAL  SHEAR  STRESS 
DISTRIBUTION 


B.  CORRECTED  SHEAR  STRESS  PROFILE 


(DEVELOPED  FROM  FIGURE  6.10) 


FIGURE  6.11  EXPERIMENTAL  SHEAR  STRESS  DISTRIBUTION 

IN  A  NON-DELAMINATED  BEAM 


SHEAR  STRESS 
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VALUES  INFLUENCED  VALUES  NOT  INFLUENCED 

BY  LOAD  PROXIMITY  ,  BY  LOAD^  PROXIMITY 


•  EXPERIMENTAL  MID-DEPTH  SHEAR  STRESS  FOR 
GAUGE  LINES  3,5,7, . ,23 


O  EXPERIMENTAL  MID-DEPTH  SHEAR  STRESS  FOR 
GAUGE  LINES  4,6,8/ . ,24 


IDENTIFIES  ’’GOOD"  RESULTS  BASED  ON 
A  CORRECTLY  SHAPED  SHEAR  STRESS 
DISTRIBUTION  AND  "X"  -  80  PSI 


-*-+  IDENTIFIES  "QUESTIONABLE”  RESULTS 
BASED  ON  A  CORRECTLY  SHAPED 
SHEAR  STRESS  DISTRIBUTION  AND 
80  PSI  -  "X"  -  200  PSI 


FIGURE  6.12 


LEGEND  FOR  FIGURE  6.13  THROUGH  FIGURE  6.20 
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6  13  MID-DEPTH  SHEAR  STRESSES  BASED  ON  ONE  GAUGE  LENGTH 

NO  DELAMINATION 
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6  14  MID-DEPTH  SHEAR  STRESSES  BASED  ON  TWO  GAUGE  LENGTHS 

NO  DELAMINATION 
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GAUGE  LINE 
(SEE  FIGURE  6.12) 


FIGURE  6  15  MID-DEPTH  SHEAR  STRESSES  BASED  ON  ONE  GAUGE  LENGTH 

DELAMINATION  LENGTH  =  0.5D 
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(SEE  FIGURE  6.12) 


FIGURE  6  16  MID-DEPTH  SHEAR  STRESSES  BASED  ON  TWO  GAUGE  LENGTHS 
h  DELAMINATION  LENGTH  =  0.5D 
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FIGURE 


GAUGE  LINE 


(SEE  FIGURE  6.12) 


6  17  MID-DEPTH  SHEAR  STRESSES  BASED  ON  ONE  GAUGE  LENGTH 

DELAMINATION  LENGTH  =  1.0D 
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SHEAR  STRESS  SHEAR  STRESS  SHEAR  STRESS 
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FIGURE 


GAUGE  LINE 
(SEE  FIGURE  6.12) 


6  18  MID-DEPTH  SHEAR  STRESSES  BASED  ON  TWO  GAUGE  LENGTHS 

DELAMINATION  LENGTH  =  l.OD 
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GAUGE  LINE 
(SEE  FIGURE  6.12) 


FIGURE  6  19  MID-DEPTH  SHEAR  STRESSES  BASED  ON  ONE  GAUGE  LENGTH 

DELAMINATION  LENGTH  =  1.5D 
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GAUGE  LINE 
(SEE  FIGURE  6.12) 


FIGURE  6.20  MID-DEPTH  SHEAR  STRESSES  BASED  ON  TWO  GAUGE  LENGTHS 

DELAMINATION  LENGTH  =  1.5D 
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FIGURE  6.21  TYPICAL  STRAIN  DISTRIBUTION  PRIOR 

TO  COMPLETE  SHEAR  FAILURE 
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FIGURE  6.22  POSITION  OF  HORIZONTAL  SHEAR  FAILURE  PLANE 
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T-18  REPRESENTS  LAMINATION  18-R 
T-20  REPRESENTS  LAMINATION  18-L 
(SEE  TABLE  6.1) 


TENSION  TEST  SPECIMENS 
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A.  PORTION  CONTAINING  DELAMINATION 


B.  PORTION  OUTSIDE  DELAMINATED  SECTION 

HORIZONTAL  SHEAR  FAILURE  SURFACE 


PLATE  6.4 


CHAPTER  VII 


DISCUSSION  OF  RESULTS 
A„  Material  Properties 
io  Moisture  Content 

TABLE  6-1  and  FIGURE  6-1A  indicate  that  the  moisture  content  of 
the  beam  decreased  and  stabilized  during  the  test  program.  The  greatest 
decrease  in  moisture  content  occurred  at  the  free  edges;  ie.  at  the  outer 
laminations  and  in  the  vicinity  of  the  delamination. 

ii.  Ultimate  Strength 

TABLE  6-1  and  FIGURE  6 - 1 B  show  that,  in  general,  the  tension 
ultimate  stress  is  greater  than  the  compression  ultimate  stress.  This 
is  particularly  true  for  the  higher  specific  gravity  values  which  occur 
in  the  outer  laminations  (FIGURE  6-2A). 

The  failure  patterns  of  the  tension  test  specimens  are  shown 
in  PLATE  6-1  and  explain  the  exceptions  to  the  above  statement. 

Lamination  specimens  9  and  16  failed  due  to  poor  grain  orientation. 

The  ultimate  strength  of  specimens  11,  14  and  19  were  reduced  as  excess¬ 
ive  crushing  at  the  grips  caused  premature  failure. 

The  ultimate  strength  characteristics  are  not  discussed  in 
further  detail  as  the  ultimate  strength  behaviour  is  not  of  prime  inter¬ 
est  or  importance  in  this  investigation. 
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iii.  Specific  Gravity 

The  specific  gravity  values  are  tabulated  in  TABLE  6-1. 

FIGURE  6-2A  shows  that  material  did  vary  from  the  highest  specific  gravity 
in  the  outer  laminations  to  a  minimum  specific  gravity  near  the  mid-depth 
lamination. 

iv.  Modulus  of  Elasticity 

TABLE  6-1  and  FIGURE  6-2B  indicate  that  the  modulus  of  elasti¬ 
city  of  Douglas  Fir  in  compression  is  slightly  higher  than  the  modulus 
of  elasticity  in  tension.  This  characteristic  has  been  noted  previously 
by  Doyle,  McBurney  and  Drow  (3).  As  the  difference  is  small  the  general 
practice  to  date  is  to  assume  equality. 

The  relation  between  modulus  of  elasticity  and  specific  gravity, 
plotted  for  each  lamination  (FIGURE  6-3),  indicates  that  as  the  specific 
gravity  increases  the  modulus  of  elasticity,  in  both  tension  and  compress¬ 
ion,  increases.  This  characteristic  has  been  reported  previously  by 
Wakefield  (14). 

A  separate  straight  line  fit,  by  the  method  of  least  squares, 
of  the  compression  modulus  values  and  of  the  tension  modulus  values 
shows  that  for  any  given  specific  gravity  the  modulus  of  elasticity  in 
compression  is  approximately  8%  higher  than  the  modulus  of  elasticity 
in  tension. 

Arithmetic  averaging  of  the  tension  and  compression  modulus 

3 

results  in  a  compression  modulus  value  of  2.05  x  10  ksi  and  a  tension 

3 

modulus  value  of  1.90  x  10  ksi.  The  average  compression  modulus  of 
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elasticity  is  therefore  8%  higher  than  the  average  tension  modulus  elasti¬ 
city. 

B.  Deflections 

Maximum  load  deflections,  based  on  a  load  1.33  times  the  working 
load,  are  recorded  in  TABLES  6-2,  6-3  and  6-4.  The  beam  remained  within 
the  elastic  range  for  all  tests.  A  comparison  of  the  maximum  load 
deflections  indicate  that  the  delamination  did  not  reduce  the  stiffness 
of  the  beam  in  the  range  of  loading  applied. 

C.  Load  Influence 

Preliminary  examination  of  the  flexural  strain  profiles  indicated 
that  the  concentrated  loads  had  significant  localized  effect  on  the  long¬ 
itudinal  strain  distributions.  As  a  result  an  examination  was  conducted 
to  determine  which  strain  profiles  were  not  effected  by  the  load  proxim¬ 
ity. 

It  was  found,  by  plotting  the  strain  profiles,  that  at  a  load 
1.33  times  the  working  load  only  gauge  lines  21  through  24  were  uneffect¬ 
ed  by  the  single  point  mid-span  load.  This  meant  that  only  four  strain 
profiles,  hence  two  shear  profiles,  were  reliable  for  the  conditions  of 
no  delamination  and  a  delamination  length  of  1.5D.  No  non-influenced 
results  could  be  obtained  for  a  delamination  length  of  0.5D  and  1.0D 
because  of  the  reduced  number  of  gauge  lines  read. 

Strain  profiles  at  working  design  mid-span  loads  were  then 
examined.  FIGURE  6-4  shows  typical  strain  profiles  for  a  span/depth  ratio 
of  8  and  no  delamination.  These  profiles  indicate  that  gauge  lines  11 
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through  24  were  not  effected  by  the  proximity  of  the  load.  The  load 
proximity  effect  was  not  changed  for  span/depth  ratios  of  10  and  12, 

The  working  load  did  not  effect  the  longitudinal  strain  profiles 
at  gauge  lines  located  at  a  distance  of  0,80  times  the  beam  depth  away 
from  the  center  line  of  the  applied  load.  When  the  load  level  was 
increased  to  1,33  times  the  working  load  the  load  effect  was  evident  up 
to  a  distance  equal  to  1.40  times  the  beam  depth. 

Similar  plots  were  made  for  the  two  point  load  conditions  (FIGURE 
6-4),  It  was  found  that  two  point  loads  effected  the  strain  profiles  in  the 
same  way  as  did  the  single  point  load. 

From  the  study  of  load  influence  it  was  concluded  that  for  the  case 
of  mid-span  working  load  the  strain  profiles  of  gauge  lines  11  through  24 
are  uneffected  by  the  load  for  the  case  of  no  delamination  and  for  the 
case  of  a  delamination  length  equal  to  1.5D.  Gauge  lines  11  through  16  and 
gauge  lines  11  through  20  are  free  of  load  influence  for  the  case  of 
delamination  length  of  0.5D  and  1.0D  respecti vely .  For  the  case  of  two 
point  loads,  at  the  working  design  load,  gauge  lines  3  through  16  were 
not  effected  by  the  loads  for  all  cases  of  span/depth  ratios. 

D.  Determination  of  Effective  Modulus  of  Elasticity 

i.  Introduction 

To  enable  a  proper  shear  stress  distribution  calculation  the 
best  possible  flexural  stress  profiles  had  to  be  defined.  A  well  defined 
flexural  stress  profile  depends  on  the  reliability  of  the  strain  data  and 
on  the  modulus  of  elasticity. 
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To  check  the  validity  of  the  test  data  the  flexural  strains  were 
examined  in  the  pure  moment  region  of  the  two  point  load  tests .  TABLE  6-5, 
typical  of  all  two  point  load  tests,  shows  the  results  of  the  best  fit 
curve  through  the  working  load  strain  values  for  the  case  of  no  delamination 
and  a  span/depth  ratio  of  8.  The  similarity  of  the  results  can  be  deter¬ 
mined  by  examining  the  two  extreme  strain  values  (TABLE  6-5).  It  should  be 
remembered,  when  examining  this  table,  that  the  strains  are  based  on  data 
read  to  an  accuracy  of  ±  0,00001 

TABLE  6-5  indicates  that  the  maximum  compressive  strains  are  about 
10%  less  than  the  tension  strains.  This  was  also  noted  by  Collins  (2), 

The  difference  in  strains  point  out  the  difference  between  the  compression 
modulus  of  elasticity  and  the  tension  modulus  of  elasticity. 

To  determine  an  effective  modulus  of  elasticity  for  the  beam  the 
fitted  strain  values  for  the  pure  moment  region  were  used.  As  the  external 
loads  and  their  location  were  known  exactly  a  correct  evaluation  of  the 
external  applied  moment  was  possible,  A  static  equilibrium  check,  based 
on  the  internal /external  moment  ratio,  was  then  made  for  each  gauge  line. 
Also  the  compression/tension  force  ratio  was  determined  as  a  second  check, 
TABLE  6-6  summarizes  the  internal /external  moment  ratios  and  the 
compression/tension  force  ratios  for  the  different  modulus  of  elasticity 
evaluated  in  the  equilibrium  check. 

ii.  E  TEST 

The  first  modulus  of  elasticity  to  be  checked  in  the  static 
equilibrium  check  was  the  modulus  of  elasticity  for  individual  laminations 
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(TABLE  6-1).  TABLE  6-6  shows  that  the  moment  ratio  check  agrees  within 
a  maximum  error  of  5%.  The  horizontal  force  check  is  not  as  good  as  the 
moment  check  particularly  for  the  tests  with  the  span/depth  ratio  of  12. 

iii.  E  BEND 

An  equilibrium  check  was  conducted  using  the  modulus  of  elasticity 
determined  from  the  mid-span  deflections.  On  the  basis  of  several  deflect¬ 
ions  the  modulus  of  elasticity  in  bending  was  found  to  be  1.73  x  10^  ksi. 

The  equilibrium  results  (TABLE  6-6)  indicate  that  the  internal 
moments  were  approximately  20%  low.  Also  the  horizontal  force  check  had  a 
larger  variance  than  the  previously  discussed  modulus  check. 

iv.  E  AVERAGE 

An  equilibrium  check  was  made  using  the  modulus  of  elasticity 
determined  by  averaging  the  modulus  of  elasticity  values  for  individual 

laminations  in  compression  and  in  tension.  This  averaging  resulted  in  a 

3  3 

compression  modulus  of  2.05  x  10  ksi  and  a  tension  modulus  of  1.90  x  10 

ksi. 

The  equilibrium  results  (TABLE  6-6)  show  an  improvement  in  the 
compression/tension  ratio  check.  A  consistent  9%  error  occurred  in  the 
check  of  the  internal /external  moment  ratio. 

v.  E  MODIFIED 

An  increase  of  approximately  10%  to  both  the  average  compression 
and  the  average  tension  modulus  values  essentially  eliminated  all  internal/ 
external  moment  errors  (TABLE  6-6).  This  10%  increase  resulted  in  a 

3  3 

compression  modulus  of  2.25  x  10  ksi  and  a  tension  modulus  of  2.08  x  10 
ksi . 
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On  the  basis  of  this  and  all  previous  equilibrium  comparisons  it 
was  decided  that  all  further  calculations  would  be  based  on  the  modified 
modulus  of  elasticity. 

E.  Shear  Stresses 

i.  Introduction 

In  a  homogeneous,  isotropic,  rectangular  beam  the  shear  stress 
distribution  across  the  depth  of  the  beam  is  a  second  degree  parabola 
having  a  maximum  value  at  the  neutral  axis.  The  average  theoretical  shear 
stress  isa  function  of  the  transverse  shear  force  and  the  cross-sectional 
area.  In  a  rectangular  beam  the  maximum  shear  stress  is  1.5  times  the 
average  shear  stress  (v  =  1.5  V/A). 

In  a  delaminated  beam  a  shear  stress  concentration  exists  in  the 
vicinity  of  the  delamination.  Therefore  the  delamination  in  the  beam  will 
cause  an  increase  in  the  shear  stress,  in  the  vicinity  of  the  delamination, 
above  the  1.5  times  the  average  shear  stress  value. 

In  the  following  discussions  the  large  variance  of  the  results 
makes  it  difficult  to  make  absolute  statements. 

ii.  A  Specific  Shear  Stress  Calculation 

FIGURES  6-10  and  6-11  show  the  development  of  a  shear  stress 
distribution  calculation  for  a  beam  with  no  delamination.  This  particular 
shear  stress  distribution  has  been  developed  between  gauge  lines  16  and  20 
for  the  case  of  mid-span  working  load  with  a  span/depth  ratio  of  10. 

FIGURE  6-10A  shows  the  distribution  of  the  longitudinal  strains 
across  the  depth  of  the  beam.  FIGURE  6-10B  indicates  the  shape  of  the  long¬ 
itudinal  stress  profile  after  the  strains  have  been  fitted  to  a  second  degree 


nofJaubotinl  . r 

fsno id 992-220*0  9f1d  bm,  99*10^  *i69il2  92n9V2n6*id  odd  "to  nordom  6  '  zzs^iz 


9fi^  29mrd  2.1  2r  229nd2  'i69ri2  murn'rxfim  9rid  m69d  "i6t upf!6:09T  6  n 


. 

grid  nr  i arx9  norlfi  tnsjnc  >  229’.d2  nfiBfla  *  rnsad  >e;  6i  rm  l9:>  6  nl 
f  r  rw  mfesd  9f1d  nr  noid6nrm&l9b  9rid  9no^9n9i1T  .nordsnrm  eb  9 4  vo  \^inro:v 
enor36nrm6f9b  srld  tfrnrorv  9r:d  nr  t229ndz  n&srtz  9dd  nr  269nonr  n*  92U&: 


229"fd2  *t6'3il2  6  to  J'P  )l  3V‘  D  9f1l  WO'4 2  *C'  .16  -12  .  .^JOII 


'i6ruD’rl'i6q  arrlT 


99^r  )  b-  )D92  6  oi  b9  Tt  roc.  r  sri  2n;6*id2  9t  ;  oils  eMlo  i  as  a  isnrbi  r 


71 


parabola  and  the  modified  modulus  of  elasticity  has  been  applied. 

FIGURE  6-11A  shows  the  first  step  in  the  mid-depth  shear  stress 
calculation.  By  taking  the  stress  difference  at  each  lamination,  between 
gauge  lines  16  and  20,  and  adding  the  individual  differences  from  the  top 
of  the  beam  to  the  mid-depth  of  the  beam  and  from  the  bottom  of  the  beam 
to  the  mid-depth  of  the  beam,  the  non-corrected  shear  stress  profile  was 
developed. 

The  shear  stress  profile  was  then  corrected  to  give  compatible 
results  at  the  mid-depth  of  the  beam.  The  corrected  mid-depth  shear 
stress  value  was  determined  by  averaging  the  two  mid-depth  shear  stress 
values  of  the  non-corrected  distribution.  The  remainder  of  the  distri¬ 
bution  was  then  corrected  by  making  proportional  adjustments,  dependent 
on  the  mid-depth  adjustment  and  the  location  of  the  point  being  corrected 
with  respect  to  the  half  depth  of  the  beam,  to  the  non-corrected  shear 
stress  distribution. 

The  resultant  corrected  shear  stress  distribution  is  shown  in 
FIGURE  6-11B.  The  "bulge"  in  the  lower  section  of  the  corrected  experi¬ 
mental  shear  distribution  profile  is  a  direct  result  of  the  correction 
method  used.  It  is  seen  that  the  corrected  shear  stress  distribution 
is  very  similar  in  shape  and  in  magnitude  to  the  theoretical  shear 
stress  distribution. 

iii.  Experimental  Mid-Depth  Shear  Stresses 

A  close  examination  of  FIGURE  6-12  will  help  the  reader  under¬ 
stand  the  terminology  and  the  reliability  of  the  results  discussed  in  this 
section.  - 

FIGURES  6-13  and  6-14  are  plots  of  the  theoretical  and  experimental 
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shear  stress  values  at  the  mid-depth  of  the  beam  for  the  case  of  no 
delamination.  Shear  values  "good"  and  "questionable"  appear,  for  the  most 
part,  to  be  well  within  an  acceptable  range  of  the  theoretical  maximum  shear 
stress  values.  As  might  be  expected  the  "good"  experimental  shear  stresses 
have  less  range  of  error  when  calculated  over  two  gauge  lengths. 

In  the  series  of  tests  conducted  for  the  delamination  length  of 
0.5D  (FIGURES  6-15  and  6-16)  the  shear  stress  values  not  effected  by  load 
proximity  fall  outside  the  range  of  the  theoretical  shear  stress  concent¬ 
rations.  The  "good"  and  "questionable"  shear  stress  values  represent 
closely  the  theoretical  maximum  shear  stress  value  based  on  no  delamination. 

FIGURES  6-17  and  6-18  are  plots  of  the  experimental  and  theoretical 
shear  stress  values  for  delamination  length  equal  to  1.0D.  "Good"  and 
"questionable"  values,  not  effected  by  load  proximity,  agree  closely  with 
Guyon's  and  Lee's  theoretical  predictions. 

In  the  case  of  the  delamination  length  equal  to  1.5D  (FIGURES  6-19 
and  6-20)  all  values  of  mid-depth  shear  stress  were  independent  of  load 
influence.  Unfortunately  very  few  points  were  classified  as  "good"  or 
"questionable".  As  result  a  comparison  of  the  experimental  and  the  theor¬ 
etical  shear  stresses  within  a  distance  of  D/4  from  the  end  of  the  delami¬ 
nation,  a  range  in  which  Lee's  and  Guyon's  theoretical  predictions  of 
maximum  shear  stress  became  widely  divergent,  is  meaningless.  However,  it 
would  appear  that  both  Lee's  and  Guyon's  theoretical  shears  are  tolerably 
accurate  a  distance  equal  to  and  greater  than  D/4  from  the  end  of  the 


delamination. 
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F.  Further  Comments 
i.  Introduction 

In  this  investigation  several  observations  require  comment.  These 
observations  include;  uniformity  of  strain  distributions,  "splitting" 
effect  of  a  delamination,  effect  of  combined  shear  and  bending  moment  on 
equilibrium  conditions,  and  the  variance  of  the  compression/tension  ratio, 
ii.  Uniformity  of  Strain  Distribution 

Collins  (2)  found  that  strain  distribution  appeared  to  satisfy 
the  normal  assumption  of  a  linear  distribution.  Ramos  (10)  found  that  the 
strain  distribution  did  contain  abrupt  variations  in  strain  through  the 
depth  of  the  section. 

In  this  investigation  the  strain  distribution  did  have  an  overall 
uniformity.  For  purposes  of  deriving  a  best  fit  curve  the  distribution  was 
assumed  to  be  a  second  degree  parabola.  A  higher  degree  curve  fit  would 
not  result  in  a  more  precise  strain  as  the  coefficient  of  the  second  power 
term  approaches  an  insignificant  value  (TABLE  6-5). 

In  general  the  strain  distribution  was  much  more  uniform  in  a 
region  of  pure  moment  than  in  regions  of  combined  moment  and  shear  (FIGURE 
6-6  through  FIGURE  6-9). 

No  strain  discontinuity  occurred  at  the  end  of  a  delamination 
under  the  pure  moment  case  (FIGURE  6-9A).  Under  the  condition  of  moment 
and  shear^  strain  discontinuity  did  take  place  in  the  vicinity  of  the 
delamination  (FIGURE  6-9B).  This  discontinuity  was  predominant  over  a 
length  of  one  gauge  interval  and  not  noticeable  at  a  distance  of  two 
gauge  intervals  (or  D/4)  away  from  the  end  of  the  delamination.  The 
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relatively  high  discontinuity  indicates  that  a  high  localized  shear 
stress,  such  as  that  predicated  by  Lee,  does  exist  near  the  delamination. 
Unfortunately  the  actual  magnitude  of  this  stress  cannot  be  quantitatively 
determined  as  pointed  out  previously  in  this  chapter  (  Section  E.  iii). 

iii.  "Splitting"  Effect  of  Delamination 

Huggins,  Aplin  and  Palmer  (5)  and  Strong  (12)  have  recorded  opposi 
opinions  as  to  how  the  delamination  behaves  as  the  load  is  increased. 

In  this  investigation  it  was  observed  that  under  the  influence  of 
single  point  mid-span  load  the  mid-span  delamination  tended  to  close. 
However,  under  pure  moment  loading  the  delamination  tended  to  open  up. 

iv.  Effect  of  Combined  Shear  and  Moment  on  Equilibrium  Conditions 

TABLE  6-7  compares  the  equilibrium  conditions  for  a  region  of  pure 
moment  and  a  region  of  combined  moment  and  shear.  The  effect  of  shear  has 
little  effect  on  the  internal /external  moment  ratio  but  produces  a  large 
variance  in  the  horizontal  force  ratio.  This  divergence  of  results  points 
out  the  anisotropic  characteristic  of  wood  and  the  complexity  of  the 
localized  stress  distribution. 

v.  Divergence  of  C/T  Ratios 

TABLE  6-6  disci oses, for  the  modified  modulus  of  elasticity,  that 
most  C/T  ratios  range  about  1.0.  Two  cases  that  diverge  greatly,  under 
pure  moment,  are  the  case  of  no  delamination  and  that  for  a  delamination 
length  of  1.5D  when  the  span/depth  ratio  was  12. 

A  close  visual  check  of  the  strain  distribution  for  the  1.5D 
delamination  case  reveals  a  distinct  upward  movement  of  the  neutral  axis 
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(FIGURE  6-9A).  A  possible  explanation  of  this  change  in  neutral  axis 
position  is  that  the  "splitting"  effect  reduces  the  curvature  in  the  top 
portion  of  the  beam,  relative  to  the  bottom  portion,  and  thereby  places  a 
reduced  loading  effect  in  the  compression  zone.  The  C/T  ratio  is  much 
improved  for  the  case  of  the  single  mid-span  load  where  separation  is 
prevented  because  of  the  loading  (TABLE  6-7). 

It  is  difficult  to  attempt  an  explanation  for  the  case  of  no 
delamination  in  which  the  compression  force  predominates  with  respect  to 
the  tension  force.  It  was  noticed  however,  for  this  case,  that  the  neutral 
axis  did  shift  down  slightly  with  no  apparent  reduction  in  the  top  lamination 
compression  stress  (FIGURE  6-6A). 

G.  Beam  Failure 

The  failure  of  the  beam  occurred  as  a  result  of  horizontal  shear. 

The  conditions  at  failure  were:  a  1.5D  length  delamination  placed  at  the 
quarter  point  of  the  span,  a  span/depth  ratio  of  8  and  a  mid-span  load. 

The  load  was  being  increased  from  no  load  to  the  working  design 
load  of  30.9  kips.  At  a  load  of  25.6  kips  it  was  observed  that  the  load 
rate  was  extremely  slow.  Examination  of  the  end  of  the  delamination 
revealed  that  a  shear  failure  had  occurred.  The  load  was  immediately 
released. 

Vertical  marks  were  placed  at  mid-depth  on  the  beam  on  either  side 
of  gauge  line  9  and  gauge  line  10.  PLATE  6-3  shows  the  subsequent  relative 
movements  between  the  top  and  bottom  portions  of  the  beam,  in  the  vicinity 
of  gauge  line  9,  for  different  load  increments.  Similar  movement  was 

observed  about  gauge  line  10. 
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At  a  load  of  15.45  kips  Demec  readings  were  taken.  A  typical 
strain  profile  (FIGURE  6-21)  indicates  that  the  beam,  in  the  instrumented 
section,  had  failed  and  had  become  essentially  two  beams. 

FIGURE  6-22  shows  the  progression  of  the  horizontal  shear  failure 
at  different  load  increments.  With  each  5  kips  load  increase  the  shear 
failure  progressed  approximately  2  feet. 

PLATE  6-4  shows  the  top  and  bottom  halves  of  the  beam  after  failure 
was  complete.  The  portion  of  the  beam  containing  the  delamination  is  shown 
on  PLATE  6-4A.  PLATE  6-4B  is  a  continuation  of  PLATE  6-4A. 
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CHAPTER  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

The  most  significant  conclusions  based  on  the  results  of  this 
investigation  are  summarized  as  follows: 

1.  Within  the  elastic  loading  range  short  length  delaminations  do 
not  reduce  the  beam  stiffness . 

2.  Load  proximity  has  a  marked  effect  on  longitudinal  strain 
distribution. 

3.  Experimental  strain  profiles  across  the  depth  of  the  beam  can  be 
satisfactorily  defined  by  using  a  Demec  gauge. 

4.  The  beam's  effective  modulus  of  elasticity  is  dependent  on  the 
elastic  properties  of  the  laminations  and  to  some  extent  on  the 
loading  condition. 

5.  Reasonably  accurate  shear  stress  distributions  can  be  computed 
from  equilibrium  considerations  using  longitudinal  strain 
distributions . 

6.  Lee's  and  Guyon's  theoretical  mid-depth  shear  stresses  are 
tolerably  accurate  at  any  section  a  distance  greater  than  one 
quarter  the  beam  depth  from  the  end  of  a  mid-depth  delamination. 
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B.  Recommendations 

Further  tests  should  be  conducted  to  determine  the  experimental 
shear  stresses  at  any  section  less  than  one  quarter  the  beam  depth  from 
the  end  of  a  mid-depth  delamination.  To  obtain  reliable  experimental  shear 
stresses  in  this  region  the  following  are  recommended: 

1.  The  load  point  should  be  at  least  the  depth  of  the  beam  from  the 
end  of  the  delamination. 

2.  A  longer  gauge  length  should  be  used. 

3.  To  obtain  shear  stresses  over  a  small  length  of  the  beam  the 
longer  gauge  lengths  should  overlap. 
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APPENDIX  A 


COMPRESSION  AND  TENSION  TESTS 

A-l  Compression  Tests 

The  cross-sectional  dimensions  of  each  test  specimens  were 
recorded  to  the  closest  0.001  inch  at  three  different  locations. 
Longitudinal  strains  were  recorded  for  each  1000  pound  load  increment. 
This  data  was  then  fed  into  the  IBM  7040  Digital  Computer. 

Final  results  included  cross-sectional  area,  values  of  stress 
and  strain  for  each  load  increment,  ultimate  stress  and  strain,  and 
the  modulus  of  elasticity  in  compression  for  each  test  specimen. 

A-2  Tension  Tests 

The  same  procedure  was  used  in  the  determination  of  the  tension 
modulus  as  for  the  compression  modulus.  Load  increments  of  500  pounds 
were  used  as  the  cross-sectional  area  was  much  less  than  that  of  the 
compression  specimens. 

Final  results  did  not  include  an  ultimate  strain  value  as 
initial  splintering  of  the  specimen  broke  the  bond  between  the  Demec 
points  and  the  wood. 
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APPENDIX  B 

DETERMINATION  OF  MOISTURE  CONTENT  AND  SPECIFIC  GRAVITY 
B-l  Moisture  Content 

Specimens  were  weighed  prior  to  being  dried  in  a  oven  at  100° 

C  ±  1°C  for  approximately  34  hours.  Prior  to  final  weighing  the  spec¬ 
imens  were  placed  in  a  desiccation  cabinet  for  2  hours.  Moisture  content 
was  then  computed. 

Moisture  content  (%)  =  (initial  weight)  -  (oven-dry  weight)  x  100 

(oven-dry  weight) 


B-2  Specific  Gravity 

Specific  gravity  was  determined  using  the  same  sample  used  in 
the  moisture  content  test.  A  container  was  exactly  filled  with  mercury 
by  pressing  a  levelling  plate  against  the  top  of  the  container.  The 
oven-dry  sample  was  then  pressed  into  the  mercury  using  the  same  level¬ 
ing  plate.  The  displaced  mercury  was  weighed  then  the  specific  gravity 
was  calculated. 

Volume  of  block  =  weight  of  displaced  mercury 

specific  gravity  of  mercury 

Specific  gravity  of  wood  =  oven-dry  weight 

volume  of  block 
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APPENDIX  C 


COMPUTER  PROGRAM  FOR  CALCULATION  OF  SHEAR  STRESSES 

The  program,  to  develop  the  shear  stress  distribution  across 
the  depth  of  a  glulam  beam  using  Demec  gauge  readings,  was  written  in 
Fortran  IV  for  the  IBM  7040  Computer. 

In  total  69,600  Demec  gauge  readings  were  made  during  the  test 
program.  As  a  result  care  had  to  be  taken  in  organizing  the  data  to 
enable  immediate  identification  of  the  data  and  of  the  results.  There¬ 
fore  all  data  cards  were  pre-identified  with  a  code.  At  the  top  of 
each  page  of  output  the  delamination  length,  delamination  position, 
load  condition  and  the  span/depth  ratio  were  recorded. 

The  program  was  further  complicated  by  lack  of  memory  space 
in  the  machine.  To  avoid  using  tapes  several  arrays  were  repeatedly 
redefined  throughout  the  program. 

The  program  was  essentially  divided  into  two  parts.  Part  I 
controlled  the  format  of  the  output.  All  computations,  with  minor 
exceptions,  were  done  in  subroutines  which  formed  Part  II  of  the 
program. 

Parts  I  and  II  are  shown  on  the  following  pages.  A  brief 
description  precedes  each  section  of  the  program. 
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PART  I  MAIN  PROGRAM 


MAIN  PROGRAM  WRITES  THE  FOLLOWING. 


A.  BASED  ON  EXPERIMENTAL  DEMEC  READINGS 

1.  DATA 

2.  MODULUS  OF  ELASTICITY 

3.  STRAIN  VALUES 

4.  PLOTS  STRAIN  VALUES 

5.  STRESS  VALUES 

6.  PLOTS  STRESS  VALUES 


B.  BASED  ON  CURVE  FITTED  STRAIN  VALUES 

1.  STRAIN  VALUES 

2.  PLOTS  STRAIN  VALUES  _ 

3.  STRESS  VALUES 

4.  PLOTS  STRESS  VALUES 

5.  EQUILIBRIUM  CONDITIONS 

6 .  AVERAGE  OF  EQUIL I BR I UM  COND I T I ON  S 

7.  SHEAR  STRESSFS  BASED  ON  1  GAUGE  LENGTH 
AND  PLOTS 

8.  CORRECTED  SHEAR  STRESSES  BASED  ON  1  GAUGE 
LENGTH  AND  PLOTS 

9.  SHEAR  STRESSES  BASED  ON  2  GAUGE  LENGTHS 
AND  PLOTS 

10.  CORRECTED  SHEAR  STRESSES  RASED  ON  2  GAUGE 
LENGTHS  AND  PLOTS 
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$JOB  180033  CARL  WILSON 

$T  IME _  8 

SIB JOB  NO GO 

SIBFTC  DECKA  DECK 

C  CALCULATION  OF  SHEAR  STRESSES  IN  A  GLULAM  BEAM.  TO  SAVE  STORAGE 

C  SPACE  SUBROUTINES  WFRF  USED.  AS  A  RESULT  THE  ST ( NGL * NLAM , NL ) 

C  ARRAY  HAS  BEEN  REDEFINED  SEVERAL  TIMES  THROUGHOUT  THE  PROGRAM. 

C  OTHER  ARRAYS  HAVE  BEEN  REDEFINED  IN  THF  PROGRAM  WHEN  REQUIRED. 

C  IT  SHOULD  BE  NOTED  THAT  EACH  SECTION,  EXCEPT  THE  DATA  READING 

C  SECTION,  Is  AN  INDEPENDENT  UNIT. 

C  8  MINUTES  ARE  REQUIRED  TO  TRANSLATE  THIS  PROGRAM,  AS  IT  APPEARS 

C  BELOW,  TO  BINARY  FORM.  THE  FOLLOWING  PROGRAM  PLOTS  CURVES  FOR 

C  ONLY  THE  WORKING  LOAD  INCREMENT.  AN  ADDITIONAL  NL  DO  LOOP  AT 

C  THE  BEGINNING  OF  EACH  PLOT  SECTION  IS  REQUIRED  IF  ALL  LOAD 

C  INCREMENTS  PLOTS  ARE  DESIRED.  ALL  VALUES  OF  ALL  INCREMENTS  ARE 

C  RECORDED  WITHOUT  THIS  LATTER  NL  REQUIREMENT.  WITH  A  BINARY 

C  DECK  10  MINUTES  AND  26000  LINES  OF  OUTPUT  ARE  REQUIRED  FOR  EACH 

C  TEST. 

C  DATA  CARDS  MUST  BE  ARRANGED  AS  FOLLOWS  FOR  EACH  TEST* 

C  IN  FRONT  OF  ALL  DATA  CARDS. 

C  1 •  COMPRESSION  MODULUS 

C  2*  TENSION  MODULUS 

C _  IN  FRONT  OF  EACH  GROUP  OF  TEST  DATA. 

C  3.  CONTROL  OF  LIMIT  AND  LOLIM 

C  4.  CONTROL  OF  LAML ,  LSPAN,  AND  LOAD 

C  5.  DEMEC  READINGS 

C 

c 

COMMON  DR (30*20, 8), ST (30, 20, 5) ,DEL ( 24 , 20 ,3 ) ,DELT ( 24 , 20 ,5 > , 

1H0RC (30,3) ,HORT(30,3>  ,SH0R(30,3>  ,RHOR(30,3) , TF I  MO ( 3 0 , 3 )  ,F T MO ( 2 0 )  , 
1 FXMO (30,3) ,SMOM(30,3> ,RMOM(30,3) ,EC ( 20 > , ET ( 20 > *P ( 20 ) , A ( 20 ) ,Y ( 20 > , 
IX ( 30 ) ,BR 1 (30,3 ) , BB2 (30,3) ,BB3( 30,3 )  , Cl (30,3  > ,C2 ( 30,3 > ,C3(30,3) , 

1 AAl  , AA2 , AA3 , AA4 , AA5 , AP2  * AP3  » AP4  * AP5  *BP 1 »BP2  *BP3 , 

1SH2,SH3,SH4,SH5,SM2,SM3,SM4,SM5,TH1 , TH2 , TH3 , TH4 , TH5 , TH6 , TH7 , TM 1 , 
1TM2,TM3,TM4,TM5,TM6,TM7,DIV,B1*B2*B3,DET, 

1B*D,CAL,TL ,GAGL , LIMIT, LOLIM, LAML,LSPAN, LOAD, PLOAD, SPAN, NGL,NLAM, 
1NL,NNL,NGLM,NGLG,NGLL ,NGLR,R ,1 , IN ,NNN, NORuN 

c 

C*******  THIS  SECTION  READS  DATA  AND  ALSO  WRITES  THE  DATA.  IT  MUST  BE 


C 

c 


- 


INCLUDED. 


B  =  3  •  0 
D=30.875 

CAL=0. 0000199 

TL=D/19 .0 
GAGL=3.90 

RFAD ( 3 , 1 000 ) (EC (NLAM)  ,NLAM  =  1 *20) 
RF AD (3, 1000) (ET(NLAM)  ,NLAM=1,20) 

3  CONTINUE 

WR I TE ( 6 , 901 0 ) 

4  READ (5 , 1001 ) LIMIT ,LOL IM 

5  REA-D(3, 1002  )LAML  ,LSPAN, LOAD 

IF ( LAML • EQ. 1 )NGLM=30 
IF (LAML.EQ.2)NGLM=22 
IF (LAML.EQ.3 ) NGLM=26 
IF (LAML. EQ. 4) NGLM= 30 
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I F ( LOAD .EQ.2)NNL=5 
IF (LOAD.EQ.3)NNL=2 
IF (LOAD.EQ.4)NNL=2 
IF(LOAD.EQ.5)NNL=2 

IF (LOLI M.EQ.4) NNL  =  4 
IF (LSPAN.EQ.l ) P LOAD =4 1.2 
IF (LSPAN.EQ.3) PLOAD=33.0 
IF (LSPAN.EQ.5) PLOAD=27. 5 
IF (LSPAN.EQ.l )SPAN=8.0*D 
I F ( L  SP  AN.EQ.3 ) SPAN  =10.0*D 

IF ( LSPAN.EQ.5 )SPAN=1 2 .0*D 

NGLG=NGLM-*6 

NGLL=NGLM-7 

NGLR=NGLM-6 

DO  8  NL=1,NNL 

DO  6  N6L*1,NGLM 

READ (5  *1004) (DR (NGL  ,NLAM,NL )  ,NLAM=1 *10 ) ,NC1 ,NC2 ,NC3  ,NC4,NC5  ,NC6 
WRITE (6, 2003) ( PR ( NGL , NLAM , NL ) ,NLaM=1 ,10) ,NC1 ,NC2 ,NC3 ,NC4,NC5 ,NC6 
RE  AD ( 5  *  1 004 )  (DR(NGL,NLAM,NL) *NLAM= 11,20) ,NC1 , NC 2 , NO ,NC4 * N CS , NC6 
WRITE  (6, 2003)  (  DR  (  NGL  ,  NL  AM  ,  NL  )  ,  NL  A.M  =  1 1 ,20)  ,  NCI  ,NC2,NC3,NC4  ,NC5,NC6 
6  CONTINUE 

WRITE  (6, 901  OF 

8  CONTINUE 

WR  ITE ( 6 , 9003 )  _  _ 

WR  I TE ( 6  *  2080 ) 

WRITE (6, 9003) 

IF(LAML.FQ.l ) WR I TE ( 6  *  2000 ) 

IF (LAML.FQ.2 ) W  R I T F ( 6 , 2002) 

I F ( L AML • EQ. 3 ) WR ITE (6,2004) 

IF ( L AML .EQ. 4) WRITE (6, 2006) 

IE (LSPAN.EQ.l)WRITE (6,2008 ) 

IF (LSPAN.EQ.3)WRITE(6 ,2012 ) 

IF (LSPAN.EQ.5) WRITE (6, 201 6) 

I E ( LOAD. EQ . 1 ) WR I TE (6,2018) 

IE (L0AD.EQ.2 ) WRITE (6 ,2020) 

IF (L0AD.EQ.3 )WR ITF (6,2022 ) 

IF ( LOAD. EQ.4) WRITE (6, 2024) 

IF (L0AD.EQ.5 ) WR ITE (6, 2 026) 

WR  I  TE ( 6 , 9003 ) 

WRITE  (6, 20821“ 

WRITE(6,2083) 

WR I TE ( 6 , 9001 ) 

DO  9  NL AM= 1 ,20 

WR  I  TE (6, 2084) NLAM, EC ( NLAM) ,ET ( NLAM ) 

9  CONTINUE  _ _ _ _ 

WRITE (6,901 OT" 

C*******  THIS  SECTION  CALCULATES  AND  WRITES  STRAIN  OF  EACH  GAUGE  LENGTH. 
C  WHEN  NOR()N=  1  STRAINS  ARE  NOT  FITTED.  N0Rl)N  =  2  STRAINS  EITTED. 

C  NORUN  CONTROL  AT  END  OE  STRESS  PLOT  SECTION. 

C  SUBROUTINES  REQUIRED 

z — -  STRATA) 

C  -  FIT 

C 

NORt)N=  1 
11  CONTINUE 

CALL  STRAIN 
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DO  10  NGL  =  1 ,NGLM 
WR  I  TE  ( 6  *9003 ) 

I F ( L AML •  EQ.  1 ) WR I TE ( 6 * 2000 ) 

IF (LAML. EO. 2 ) WR I TF ( 6  * ? 002 ) 

IF (LAML.EQ.3 ) WRITE (6  *2004) 

_____  I F  (LAML  . FQ . 4 ) WR I TE ( 6  *  2 006 ) 

IF (LSPAN.EO.l )WPITF (6*2008 ) 

IF (L  SP AN .EQ. 3 ) WR ITE(6*2012) 

IF(LSPAN.EQ.5) WRITE (6*2016) 

IF ( LOAD. FQ . 1 ) WR I f F ( 6 *2 0 1 8 ) 

IF (L0AD.EQ.2 )WRITF (6  *2020) 

IF (L0AD.FQ.3 ) WR ITF (6  *2022 > 

IF (L0AD.EQ.4)WRITE(6*2024) 

IF ( LOAD. E0.5) WRITE (6 *2026) 

WR ITE ( 6  *9003 ) 

WRITE (6*2030) NGL 
WR I TE ( 6  * 0001 ) 

WR I TE ( 6  *  2032 ) 

IF (N0RUN.EQ.2 )WRI TE (6  *3004) 

WRITF(6*9003) 

I F ( NNL • GT • 2 ) GO  TO  12 
I F ( NNL •  EQ  •  2 ) GO  TO  15 
12  WR ITE(6*2034) 

WR ITE ( 6 ♦ 9001 ) 

DO  14  NL  AM  = 1*20 

WRITE (6  *2036) NLAM , (ST (NGL, NLAM, NL )  ,  NL  =  2  *  NNL ) 

14  CONTINUE 
GO  TO  18 

15  WR ITE(6  *2042 ) 

WR  I  TE ( 6  *  9001 ) 

DO  16  NLAM= 1*20 

WRITE (6 *2044) NLAM, ( ST ( NGL ,NLAM ,NL ) * NL  =  2  *NNL ) 

16  CONTINUE _  _ 

18  CONTINUE 

I F ( NORUN .EQ. 1 ) GO  TO  17 

WR  I  TF ( 6  *9003 )  _ 

WRITE (6 ♦30001 
WR I TE ( 6  *  9003 ) 

DO  17  NL-2  *NNL 

WR ITE (6  *  3002 ) NL  *C1 ( NGL  *  NL ) ,C2 ( NGL ,NL ) , C3 ( NGL , NL ) 

WRITF(6*9001 ) 

17  CONTINUE 
WRITE (6  *90101 

10  CONTINUE 
C 

Q*******  THIS  SECTION  PLOTS  THE  STRAIN  VALUE  OF  EACH  GAUGE  LENGTH. 

SUBROUTINES  REOUIRFD 
PLOT 

DO  20  NGL=1 ,NGLM 
NL  =  4 

IF (LAML.EQ.l ) WR I TE ( 6  *  2000 ) 

IF (LAML.E0.2) WRITE (6*2002 ) 

I F ( LAML • FQ. 3 ) WR ITE (6*2004) 

IF ( LAML .F0.4 ) WR ITF(6*2006) 

IF (LSPAN.EO.l )WR ITF (6 *2008 ) 

IF ( LSPAN.EQ.3 ) WRI TE (6 *2012 ) 

!P(LsPaN.E0.5  )  WRite  (6  *2oi^>  i 
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I F ( LOAD •  EQ. 1 )WRITE(6*2018) 

IF ( LOAD • FQ.  2 ) WR I^F(6*2020) 

IF (L0AD.FQ.3 ) WR ITF (6  *2022 ) 

IF ( LOAD .  FQ.  4 )WR ITF (6*2024 ) 

IF ( LOAD. EQ. 5) WRITE (6 *2026) 

_ WRITE (6 *9001  ) 

WRITE(6*2030)NGL 
WR ITE (6  *4203 ) 

IF ( NORUN .EQ .2) WRITE (6 *3006 ) 

WRITE (6 *4266) 

WRITE(6*4202) 

WR  ITE ( 6  *4204 ) 

WR I TE ( 6  *  9001  ) 

NLAM=1 
26  CONTINUE 

P  (  NLAM  )  »ST  (  NGL  , NL AM  ,NO*l 000000.0 
CALL  PLOT 

IF (NLAM.EQ.IO)GO  TO  28 
IF(NLAM.GT.10)G0  TO  30 
WR ITE ( 6  *9001 ) 

NL  AM  =  NL AM+l 
GO  TO  26 

28  CONTINUE 
NLAM=NLAM+1 
GO  TO  26 
30  CONTINUE 

WR ITE (6  *9001 ) 

NL AM=NL AM+l 

IF (NLAM.GT.20) GO  TO  32 
GO  TO  26 
32  CONTINUE 

WR ITE ( 6  *4204 ) 

WR  ITE(6*4202) 

WR ITE (6 *4 26 6 ) 

WRITE(6*2038)NGL,NL 
WRITE (6 *9010) 

22  CONTINUE 
20  CONTINUE 
C 

C*******  Tffis  SECT  TOW  CALCOLATESANd  WRITES  STRESS  OF  EACH  GAUGE  LENGTH. 


C  WHEN  NORUN= 1  STRESSFS  NOT  FITTED.  N0RUN=2  STRESSES  FITTED. 

C  SUBROUTINES  REQUIRED  _ 

C  STRFSS 

C  FORCE 

C  MOMENT 

“T  ~"STRA I N 

C  FIT 

C  STRESS 

C  ~~~ 


CALL  STRESS 
CALL  FORCE 

~TF (LOAD.GT .2 ) GO  TO  41 
CALL  MOMENT 
41  CONTINUE 
CALL  STRAIN 

IF ( NORUN .EQ. 2 ) CALL  FIT 
CALL  STRESS 

- DU  40  NGL=T*TTGTF!  ” 
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WRITE(6*9001  ) 

IF(LAML.EQ.l ) WRITE (6 *2000) 

IF (LAML. EQ.2) WRITE (6 *2002) 

I F  (LAML.EQ .3 ) WR ITE  ( 6 *2004) 

IF (LAML.EQ.4)WR ITE (6  *2006 ) 

_  IF (LSPAN.EQ.l )WRITE (6*2008 ) 

IF (LSPAN.E0.3) WRITE (6  *2012 ) 

IF(LSPAN.EQ.5)WRITE(6*2016) 

I F ( LOAD •  EQ.  1 )  WR I TE (6*2018) 

IF ( LOAD. EQ.2 ) WR  ITE (6  *202  0) 

IF ( LOAD. FQ. 3) WRITE (6  *2022 ) 

IF (L0AD.EQ.4JWR ITE (6*2024) 

IF (LOAD.EQ.5 )WR ITE (6  *2026 ) 

WR I TE ( 6  * 9001 ) 

WRITE(6*2030)NGL 
WRITE (6 *9001 ) 

WR  ITE(6*2033) 

IF (NORUN. EQ.2 ) WRITE (6  *3004 ) 

WRITE (6 *9001 ) 

IF ( NNL.GT .2 ) GO  TO  42 
IF (NNL.E0.2 )GO  TO  45 
42  WR I  TE ( 6  *  2034 ) 

WR ITE ( 6  *9001 ) 

DO  44  NLAM=1*20 

WRITE(6*2037)NLAM, ( ST ( NGL *NLAM*NL )  *  NL  =  2 ♦ NNL ) 

44  CONTINUE 
GO  TO  48 

45  WR ITE ( 6  *  2042 ) 

WR  I  TE ( 6 ♦ 9001 ) 

DO  46  NL AM= 1  *  20 

WRITE(6*2045)NLAM, ( ST ( NGL * NL AM *NL )  ,NL  =  2  *NNL) 

46  CONTINUE 
48  CONTINUE 

WR ITE ( 6 ♦ 9003) 

DO  52  NL  =  2  *NNL 
WRITE(6*2038)NGL*NL 

WRITE(6*2040)H0RC(NGL*NL) ,HORT ( NGL ,NL ) *SHOR (NGL *NL ) ,  RHOR ( NGL * NL ) 
IF(L0AD.GT.2*0R.NGL.GT.NGLG)G0  TO  54 

WRITE (6 *2046 ) TF  I  MO ( NGL *NL ) ♦ EXMO ( NGL  *NL ) , SMOM ( NGL ,NL )  , RMOM ( NGL , NL ) 
54  CONTINUE 

WRITE (6 *9001) 

52  CONTINUE  _ _ 

WR I TE (6*9010) 

40  CONTINUE 

C  _  _  _ _____  _ _ _ _ __  _ _  _ 

C*******  THIS  SECTION  PLOTS  THE  STRESS  VALUE  OF  EACH  GAUGE  LFNGTH. 

C  SUBROUTINES  REQUIRED 

C _  plot 

c 

DO  60  NGL=1,NGLM 
NL  =4 

I F ( L  AML.EQ.l ) WRITE (6*2000) 

IF(LAML.EQ.2)WRITE(6*2002) 

IF (LAML.EQ. 3 ) WR I TE ( 6  *  2004 ) 

I F  (  L  AML  .  FQ  .  4  )  WR  I  TF  (  6  *  2 0 06  ) 

IE ( LSPAN.EQ. 1 ) WRITE (6 *2098 ) 

^TF  (  LSPAN.EQ.3  )  WR  ITE(6  *2012  ) _ _ _ _ _ _ 

- IF(L5PAN.EQ.5  > WRITE (6*zoib J 
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IF ( LOAD.EQ. 1 ) WR I TE ( 6  *  2018  ) 

IF (LOAD.EQ. 2 )WR I TE (6 *2020) 

I F ( LOAD  • EQ. 3 ) WR I TE (6*2022) 

IF ( LOAD • EO . 4 ) WR I TE ( 6*2024 ) 

I F ( LOAD • FQ.  5 ) WR I TE ( 6  *  2026 ) 

WR I TE ( 6  *  9001 ) 

WR  ITE(6*2030) NGL 
WR ITE (6  *4201 ) 

IF ( NORUN .EO. 2) WRITE (6, 8006) 

WRITE (6 *4266) 

WRITE(6*4202) 

WR I T E ( 6  * 4204 ) 

WR I TF( 6*9001) 

NL AM=1 
66  CONTINUE 

P ( NL  AM) = ST ( NGL ,NLAM,NL) *0*50 

CALL  PLOT 

I F C NL AM >  EQ  > 1 0 ) GO  TO  68 
IF (NLAM.GT.IO)GO  TO  69 
WR I TE ( 6  *  9001 ) 

NLAM=NLAM+1 
GO  TO  66 

68  CONTINUE 
NLAM=NLAM+] 

GO  TO  66 

69  CONTINUE 

WR I  TE ( 6  *  9001 ) 

NLAM=NLAM+1 

IF (NLAM.GT.20)G0  TO  70 
GO  TO  66 

70  CONTINUE 

WR ITE (6  *4204 ) 

WR  I  TE ( 6  *  4202 ) 

WRITE (6 *4266) 

WRITE(6*2038)NGL *NL 

WRITF(6*2040)H0RC(NGL  *NL) *HORT(NGL,NL ) ,SHOR (NGL,NL ) ,RHOR( NGL  *NL ) 
IF(L0AD.GT.2«0R.NGL.GT.NGLG)G0  TO  79 

WRITE(6*2046)TFIM0(NGL*NL) *EXMO(NGL  *NL )  ,SMOM( NGL  *NL ) , RMOM ( NGL , NL ) 
79  WR I TE ( 6  *  9010 ) 


60  CONTINUE 

NORUN =N  OR  UN +  1  _ _ 

IF (NORUN.LT. 3 )G0  TO  11 
C 

c*******  THIS  SECTION  WRITES  EQUILIBRIUM  CONDITIONS  FOR  ENTIRE  TEST. 
C  SUBROUTINE  S  RFOUTRFD 

C  FORCE 

C  MOMENT 

CALL  FORCE 

IF ( L0AD.GT.2 ) GO  TO  81 
CALL  MOMENT 

— errrwrimrr" 

IFILAML.FQ.l ) WR I TF ( 6  *  2000 ) 

IF  ( LAML • EQ. 2 ) WR I TF  (6*2002) _ _ _ 

IF (LAML.EQ.3 )WRITF(6  *29  04) 

IF (LAML.EQ.4)WR ITE (6  *2006 ) 

IF (LSPAN.EQ.l )WRITE(6 *2008 )  _ 

- IF (LSMAN.LU.3 J WRITE (6*2H1TJ 
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IF(LSPAN.EQ.5)WRITE<6*2016> 

IF (L0AD.EQ.1 )WR ITF (6  *2018 ) 

IF (L0AD*EQ.2 )WRITF (6  *202  0) 

IF ( LOAD •FQ.3 ) WR ITF (6  *  2022 ) 

IF (L0AD#EQ*4) WRITE (6*2024) 

IF (L0AD.EQ.5 ) WR ITE(6  *2026 ) 

WR ITE(6 *9001 ) 

WR I TF ( 6  *  205  2 ) 

WR ITE( 6  *  2054 ) 

WR  I  TE ( 6  *  2056 ) 

WR ITF ( 6  * 9001 ) 

DO  80  NGL  =  1 *10 
DO  82  NL=2*NNL 

WRITF(6*2058)NGL*NL*HORC (NGL *NL) *HORT ( NGL  *NL ) , SHOR ( NGL , NL ) * RHOR ( NG 
1L,NL) *TFIMO(NGL,NL) *FXMO ( NGL ,NL ) * SMOM ( NGL *NL ) ,RMOM ( NGL * NL ) 

8?  CONTINUE 

WR I TE ( 6  * 9001 ) 

80  CONTINUE 

WR I TE  C 6*9010) 

I F ( L AML •  EQ#  1 ) WR I TE ( 6  *  2000 ) 

IF (LAML.EQ.2 )WR ITF (6  *2  002) 

I F ( L  AML • EQ . 3  )  WR I  TE  ( 6 ♦ 2004) 

IF(LAML.EQ#4)WR ITE(6*2006 ) 

IF (L5PAN.E0.I ) WRITE (6*2008) 

IF (LSPAN#EQ.3)WRITE(6*2012) 

IF(LSPAN.EQ#5)WRITE(6*2016) 

IF (LOAD.EO.l )WRITF(6* 2018 ) 

IF ( LOAD • FO • 2 ) WR ITF (6  *  2020 ) 

IF (L0AD#EQ#3 ) WR ITF (6 *2022) 

IF (LOAD* EQ. 4 ) WR ITF (6  *2024) 

I F ( LOAD #FQ#5)WRITE (6*2026) 

WR ITF ( 6  *9001 ) 

WR ITF ( 6  » 2052 )  _  _ _ _ 

WRITE (6 *2054) 

WR  ITF ( 6  *2056 ) 

WRITE (6 *9001 ) 

DO  84  NGL=1 1  * 20 
DO  86  NLS2  *NNL 

IF(L0AD.GT.2.0R.NGL.GT.NGLG)G0  TO  88 

WRlTE(6 *2058) NGL *NL*HORC (NGL *NL) *HORT ( NGL *NL > *SH0R(NGL*NL ) *  RHOK ( NG 
1L,NL) ,TFIMO(NGL,NL) » E XMO ( NGL *NL ) , SMOM ( NGL *NL ) *RMOM ( NGL * NL ) 

GO  TO  86  _ 

88  CONTINUE 

WRITE! 6 *2058) NGL  *NL *HORC (NGL *NL)  ,HORT ( NGL * NL >  * SHOR ( NGL *NL )  > RHOR ( NG 
1L  *NL ) 

86  CONTINUE 

WR I TE ( 6  *  9001 ) 

84  CONTINUE 

WR  I TF ( 6  *  9010 ) 

IF (LAML.EQ.l ) WRITE (6  *2000) 

IF (LAML.EQ.2 ) WR ITE (6  *2002  > 

IF (LAML.EO. 3 ) WR I T  F ( 6 ♦ 2  0  0  4  > 

IF (LAML#EQ.4)WR ITF (6  *  2006 ) 

IF(LSPAN.FQ#1 )WRITE(6*2008)  __ 

IF(LSPAN.EQ.3)WRITF(6  *2012 ) 

IF(LSPAN#FQ#5 ) WRITE (6*2016) 

IF(LOAD.FQ.l ) WRITE (6  *2018 ) 

TFTrOATD‘ •  TO.TTW  I  TE  (6*2020) 
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IF (L0AD.EQ.3 )WR ITF (6  *202  2 ) 

I F ( LOAD • FO • 4 ) WR ITE (6 *2024) 

IF (LOAD.EQ.5 ) WR ITF (6  *  2026 ) 

WRITFT6 *  9001 ) 

WR I TF ( 6  *  205  2  ) 

WR ITF ( 6  * 2054 ) 

WR ITF ( 6  *  2056 ) 

WR I TE ( 6 ♦ 9001 ) 

DO  90  NGL=21 *NGLM 
DO  92  NL  =  2  *NNL 

IF (L0AD.GT.2«0R.NGL.GT.NGLG)G0  TO  94 

WRITF(6*2058)NGL*NL*HORC (NGL,NL) , HORT ( NGL , NL >  * SHOR ( NGL * NL ) *  RHOR ( NG 
1L.NL) *TFIM0(NGL,NL) *EXMO(NgL ,NL ) ,SM0M( NGL.NL )  ,RM0M ( NGL  ,  NL ) 

GO  TO  92 
94  CONTINUE 

WR ITE ( 6 ♦ 205  8 ) NGL  *NL  *HORC ( NGL ♦NL ) *HORT ( NGL  *NLTT, SHOR ( NGL  *NL ) *RH0R ( NG 
1L  *NL ) 

92  CONTINUE 

WR I TF ( 6  *  900 1  ) 

90  CONTINUE 

WR I TE ( 6  *  9010 ) 

C 

c*******  THIS  SECTION  WRITES  A  SUMMATION  OF  C/T  AND  INT/EXT  RATIO 
C  VALUES.  IT  MUST  FOLLOW  SUBROUTINE  FORCE  AND  SUBROUTINE  MOMENT. 

C  IT  HAS  BEEN  INCLUDED  TO  ENABLE  A  QUICK  VISUAL  CHECK  OF  THE 

C  EFFECT  OF  THE  MODULUS  OF  ELASTICITY  ON  THE  CONDITIONS  OF 

C  EQUILIBRIUM. 

C  SUBROUTINES  R'EQU  I  RED 

C  COMPAR 

c 

I F ( NNL • EQ. 2 ) GO  TO  127 
WR ITE ( 6 ♦ 2080 ) 

WR I TE ( 6 ♦ 9001 ) 

TFTL AML .FQ. 1TWR I  TE (6  *  2000) 

I F  ( L  AML • FQ. 2 ) WR ITE(6*2002) 

IF(LAML.EQ.3)WRITE(6*2004) 

IF(LAML.FQ.4)WRITF(6*2006) 

IF (LSPAN.EQ.l ) WRITE (6 ,2008 > 

IF ( LSP AN. EQ. 3) WRITE (6*2012  > 

IF(LSPAN.EQ.5)WRtTF (6*2016 T 
IF ( LOAD .EQ. 1)WRITE(6*2018) 

IF (L0AD.EQ.2 )WR ITE (6 *2020) 

IF (L0AD.FQ.3)WRITF (6  *2022) 

IF (LOAD. EQ. 4 )WR ITF (6  *2024) 

IF(LOAD.FQ.5)WRITE(6*2026) 

WRITE (6 *9003 ) 

WR ITE(6*2082) 

WR ITE(6*2083)  _ 

WR I TF ( 6 ♦ 9001 ) 

DO  129  NLAM=i,20 

WRTTE(6*2084)NLAM,EC(NLAM) *ET(NLAM) 

— rT^xoTrnwE 

CALL  COMPAR 

WR I TE ( 6  ♦  9003 ) _ 

IF (LOAD.EQ.l ) WRITE (6  *2090)NGLG 
IF(LOAD.EQ.2)WRITF(6*2091 ) 

WR  ITE ( 6  *9003 ) 

— —  WRI  It  (6*20920 
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C-ll 

WR I TE ( 6  *  2094 ) 

WR I  TE  (  6  *9001 ) 

WR ITF(6 ,2096) THl ,TH2 ,TH3 *TH4, TH5 ,TH6 ,TH7 
WR ITE ( 6  *9001 ) 

WRITE (6 *2098 ) TM1 * TM2  * TM3  * TM4 * TM5 , TM6 * TM7 
WR  I  TE  (  6  *  901  0  ) 

127  CONTINUE 
C 

C********THIS  SECTION  CALCULATES  AND  WRITES  SHEAR  FORCES  FOR  EACH 
C  LAMINATION.  WHEN  NORUN=l  SHEARS  ARE  FOR  SUCCESSIVE  GAUGE  LINES. 

C  WHEN  NORUN=2  SHEARS  ARE  CORRECTED  AND  ARE  FOR  SUCCESSIVE  GAUGE 

C  LINES.  WHEN  NORUN«3  SHEARS  ARF  FOR  ALTERNATE  GAUGE  LINES.  WHEN 

C  N0RUN*4  SHEARS  ARE  CORRECTED  AND  ARE  FOR  ALTERNATE  GAUGE  LINES. 

C  THE  NORUN  CARRIES  THROUGH  THE  PLOT  OF  TOTAL  SHEAR  DISTRIBUTION. 

C  SUBROUTINES  REQUIRED 

C  STR A I N 

C  FIT 

C  STRESS 

C  FORCE 

C  DELTA 

C  CORREC 

C 

NORtJN  =  l 
140  CONTINUE 

CALL  STRAIN 
CALL  FIT 

_ CALL  STRESS _ _ _ _ 

CALL  FORCE 
I F ( NORUN .LE • 2 ) 1=3 
IF (NORUN.LE .2 ) IN  =  2 
IF ( NORUN.GE .3 ) 1=5 
IF (NORUN. GE. 3 ) IN=4 
CALL  DELTA 

IE (NORUN. EQ. 2 1  CALL  CORRFC 

IF ( NORUN .EQ. 4 ) CALL  CORREC 

WR I TE ( 6  * 901 0 )  _ _ _ _ _ 

DO  100  NGL= I *NGLG 

IF ( LAML .EQ. 1 ) WR I TE ( 6  *  2000 ) 

IF (LAML.EQ.2 )WR ITE (6 *2002 ) 

“TFT LAML.EQ.3 ) WR ITE (6*2004 j 

I E ( L  AML • EQ. 4 ) WR ITF(6*2006) 

IF (LSPAN.EQ.l ) WRITE (6  *2008 ) 

IF (LSPAN.EQ.3 ) WRITE (6  *2012) 

IF(LSPAN.EQ.5)WRITE(6*2016> 

IF ( LOAD.EQ. 1  )WR ITE (6  *  2018 ) 

IF ( LOAD . EQ • 2  >  WR I TE ( 6  *  2 02 0 ) 

IF ( LOAD. EQ. 3 ) WRITE (6  *2022  ) 

IF (LOAD.EQ. 4) WRITE (6 *2024) _ 

IF ( LOAD.EQ. 5 )WR ITF (6  *  2  026 ) 

WRITE (6 *9002) 

IF (NORUN. LE.2 )NNN=NGL-2 
I F ( NORUN • G  E • 3 ) N NN  =  NGL-4 
WR ITE(6*2086)NNN*NGL 
WR I TE ( 6  *  9001 ) 

WRITF(6*2031 ) 

WR ITE ( 6  *9001 ) 

IF (NNL.GT.2 )G0  TO  102 

- - - - IF  (NNL. LU#2  )L>U  !U  105"” 
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102  WR  ITE ( 6 ♦ 2034 ) 

WRITE (6, 9001 ) 

DO  104  NL AM= 1 >20 

”  WR ITE  (6 *2037)NLAM, ( DEL { NGL , NL AM ,NL )  ,NL  =  2  *NNL ) 

104  CONTINUE 
GO  TO  108 
“TO  5  WRITE(6 *2042 ) 

WR I TE ( 6  * 9001 ) 

DO  106  NLAM= 1*20 

WR  I  T  FT  6  *  2  0  4  5  )NL  A  M  ,  { DTLlTlGTr,7T[^  ,Wr,TMlT=  2  *N  NET 
106  CONTINUE 
108  CONTINUE 

WR I TF ( 6 ♦ 9008 ) 

WR I TE ( 6  ♦  2035 ) 

I F  (  NORUN  •  EQ  •  2  •  OR  •  NORUN • EQ. 4  )  WR  ITE(6*2039) 

WF I T E ( 6  *  9001) 

IF(NNL.GT.2)G0  TO  112 
I F ( NNL.  EO#  2 ) GO  TO  118 
112  CONTINUE 

DO  114  NLAM=1 *20 

WRITE(6*2037)NLAM, (DELT(NGL,NLAM,NL) ,NL=2*NNL) 

114  CONTINUE 

GO  TO  118 

115  CONTINUE 

DO  116  NLAM  =  1 *  20 

WRITE(6*2045)NLAM, ( DFLt ( NGL ,NLAM ,NL ) , NL  =  2  *NNL ) 

116  CONTINUE 
118  CONTINUE 

WRITF (6  *9010) 

100  CONTINUE 
110  CONTINUE 
C 

C********THIS  SECTION  PLOTS  THE  TOTAL  SHEAR  DISTRIBUTION.  THE 

C - DISTRIBUTION  DEPENDS  ON  NORUN  AS  EXPLAINED  TN  THE  PREVIOUS 

C  SECTION. 

C  SUBROUTINES  REQUIRED 

C  PLOT 

C 

DO  120  NGL=5 ♦ NGLG 
NL=4 

I F ( L AML • EQ. 1  ) WR I TE ( 6  *  2000 ) 

IF  (LAML.EQ.2  )WR  ITE  (6  *2002 
IF (LAML.EQ.3 ) WR I  TP ( 6  *  2004 ) 

I F ( L AML .EQ.4 ) WR I TE ( 6 ♦ 2006 ) 

IF (LSPAN.FQ.l ) WR I TE ( 6  *2008 )  _ 


IF (LSPAN.FQ.3)WRITF (6*2012) 
IF(LSPAN.EQ.5)WRITE (6  *2016) 

IF ( LOAD.EQ. 1  )WR ITF (6  *  2018 ) 

IF (LOAD. EQ. 2 ) WRITE (6  *202  0) 

IF (LOAD.EQ. 3 ) WRITF (6*2022  > 

IF (LOAD .EQ.4) WRITE (6 *2024) 

IF ( LOAD. FQ. 5) WRITE (6 *2026) 

WR I TF ( 6  *9001 ) 

IF (N0RUN.LE.2 )NNN=NGL-2 

IF ( NORUN .GE. 3 ) NNN=NGL-4 

WR  ITE(6  *2086 ) NNN *NGL 

WR ITE(6  *9001  ) 

- — —  WRITE  (T*  A20'8'")  NT” 

<tr-3 


(isr^c«A)lTI  fl’**  ^0  r 
(  foro,^)^TT  9W 
^  C  *  r ->/  .\  ji/  in  r 

•  J  .  .  J  '  .  '  -  )  ,  'A  Jl'\  (- r  r>C  •  A  )  -*>T  T  ■■ 

apt 

POf.  OT  OO 
(  'V-or  #  A  )  q-r  T  CM  3  f 
(  fono  •  X  )  IT  I 

oc,r~vMn  aor  on 

(’  .  V'=  '.('  .  I'*.lrinj^0  ♦’•A  J/(  r-  A^C,  ^  )  -IT’  r.l  1 

1'IIAIT^O')  AOt 

iumitmoo  ant 

(  fro,  A  )  ITT  cw 

( ?? ns# a) gw 

'prnC'*)  ;  j  j  qw  (  A  «0 3  •  I'M  !/K  !■< .  9f  .  S  .  03  •  Ml  90H  >  1 1 

(  r^oo#  A ) ^TT  flw 

srr  nT  on  c.to.  jhm>  qi 

arr  0T  00(S.0?#JHM)qi 

iumitmod  srr 
os .  r=MAjn  Af  r  on 

(  '  "  ■/  .  S  =  Ji\  (  _M#  MA  Jtf .  J  ■  )  T  Jtn  )  ,  ^  AJM  (trnc:,  ,>)1TT  9'-' 

qi'MTTHOO  Mr 
a  r  r  ot  oo 
iumithoo  <*  r  r 
os  •  rsMA  jia  a  r  r  on 

(  ji-v/.s.  =  jr.  (  _M.  ma  jh.  j  "v5 '  7  '  "»n ) .  VAJt*  m  a^s  •  a  )  it  1 9 

mi'/TTMoo  *rr 
9UMITM03  MT 
[Ofoo,  A)1TT  Q'.« 

^Ht/jT^oo  onr 
3UMITH03  nrr 

0 

3HT  ••  T"  fflI9T?I0  9A^HP  J  a  ^HT  ^ T  |TDB  8  2IHT***###**D 

?!U!'  :-n  h.t  ^  t  -1/-IAJ9X1  A  A  HIJ90K  tfO  ?0M9990  MOlTIJfl  ’ 9T? 10  3 


I/O  T  TAi? 


3 

3 

3 

3 


099111030  '’3MITi!OCPif? 

TO  J9 


OJdH#?  =  JrOM  ncr  on 

A=  JlA 


(  Arc,  M  IT  I  9W  f  r  .0^.  JMA  J )  9T 
Uaac,  A)-ir  9V(  r  .0“  •  J^A  J ) 9  T 

•  AOOS ♦ d)BTIflWf  t«03» JMA  I) 

(  aoot  .  a  )  tt  I  9W(  a. 09#  jma  htt 

(  'one,  ,v  )7TJ  -  v(  r.Ol.l/Aqp  J)  9T 
(  9  r0S*  A)  1TT9  •  <  J>  9? 

f  A  rr S*  A)  7T!  ^  (  9.03.M/fl ?  ))  9! 

*fOS#d>  B T I  9 w  (  r*03»CM 

(Arc^iiji  qw  (  S  •  O'3  •  0  AO  J)  q  T 
< ?  sos ,  a ) it ! fly ( r.oi.nAOj )  9T 
(  A  <*nc:  ,  A)  9T  J  p>;(  a.O^.oad  dm 
(  vror  .  -jtJ  PV.f(  P.  0^.0  AO  'I 


(  rooo,  itT  9'.* 


io,a-i/'/  i(r #  .Miiflrt/) ^ t 


A—  I  "r'/  / '/  f  P.30#!AI,rj0H) 


J0|/»  / tl M  (ARAC  .  a  1  7TT  qv 

(  f  AAf>,  A  )  OT  I  PM 
j  4]  T  ]  <»  .• 


013 


WR I TF ( 6  * 9001 ) 

WR ITE ( 6  *2035 ) 

IF (N0RUN.EQ.2.0R.N0RUN.EQ.4) WRITE <6 *2039 ) 

IF ( NOR  U  N • EQ .1 . OR . N  OR  u N . E  Q . 3 ) W  R I TE (6*9 001 ) 

WRITE(6*4?06) 

WRITF(6  *4246) 

WRITF(6  *420? ) 

WR I TE ( 6  * 4204 ) 

WR ITE ( 6  *9001 ) 

WR ITE ( 6  *4246 ) 

NLAM=1 

126  CONTINUF  _  _  _ _ 

P(NLAM) =-(DELT(NGL*NLAM*NL) ) *3.0+600.0 
CALL  PLOT 

IF ( N  L  A  M • EQ • 1 0 ) GO  TO  12R 
IF (N L AM . G T .10) GO  T 0  1  30 
WRITE (6 *9001 ) 

NL AM=NL AM+1 
GO  TO  126 
128  CONTINUF 

NL AM=NL AM+1 
GO  TO  126 
130  CONTINUE 

IF (NLAM«EQ.20)WRITF (6*4246) 

WRITE(6*9001 ) 

NL AM=NL AM+1 

IF(NLAM.GT.20)GO  TO  132 
GO  TO  126 
132  CONTINUE 

WR I TE ( 6  *4204 ) 

WRITE(6*4202) 

WRITE (6 *4246) 

WR I TE ( 6  *  901 0 ) 

122  CONTINUE 
120  CONTINUE 

NO RUN = NORUN +1 
IFCNORUN.LT. 5)G0  TO  140 
C 

C*******  THIS  SECTION  CONTAINS  THE  FORMAT  STATEMENTS  FOR  THE 
C  MAIN  PROGRAM. 

C 

IF (LIMIT. EO.O )G0  TO  3 

1000  FORMAT (10X*5F10.0*20X/10X*5F10.0*2nX/10X*5F10.0*20X/inX*5F10.0*?nx 
1  ) 

1001  FORMAT(  8X*  I2*8X*  I? )  __  _____  _ 

1002  FORM AT ( 8 X  *  12  * 8X  *  12* 8 X  *  1 2  ) 

1004  FORMAT ( 5X *  1 0F5 . 0  *  1 2X *  1 1*1X*  I2*1X*I1*1X*I1*1X*I2*1X*I1) 

2000  FORMAT ( 1X*16H  NO  DFL AM  I N AT  I  ON ) 

2002  FORMAT ( 1 X *27H  DELAMINATION  LENGTH  =  0.5D) 

2003  FORMAT  (  5X  *10F9.0*12X  *  I1*1X*I2*1X*I1*1X*U*1X*I2*1X*I1) 

2004  FORMAT ( IX  *27H  DELAMINATION  LENGTH  =  1 .OP) 

2006  FORMAT ( 1 X ♦ 2  7H  DE  L AM IN  ATI  ON  LENGTH  =  1.5D) 

2008  FOR-MAT  (  2X  *21 H  SPAN/DEPTH  RATIO  =  8) 

2012  FORMAT  (  2X  *2?H  SPAN  /DFPTH  RATIO  *  10)  __ 

2016  FORMAT ( 2X ♦ 22H  SPAN/DEPTH  RATIO  =  12> 

2018  FORMAT ( 3X * 39H  LOAD  AT  CENTRE*  DELAMINATION  AT  CENTRE) 

2020  FORM AT ( 3X * 39H  TWO  POINT  LOAD*  DELAMINATION  AT  CENTRE) 
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2024  FORMAT ( 3X  ,  39H  LOAD  AT  CENTRE  ,  1/4  POINT  DELAMINATION) 

2026  FORMAT ( 3X,42H  LOAD  AT  3/4  POINT*  1/4  POINT  DEL AM  I N A T I ON ) 

2030  FORMAT ( 7X*11H  GAUGE  LINE*I3) 

2  031  FORM  AT  (  3  OX  *  2  OH  SHF  AR.  STRFSS  I N  PS  I  ) 

2032  FORMAT ( 2  9X  ,  2  3H  STRAINS  IN  INCHES/INCH) 

2033  FORMATt 28X,24H  STRESSES  IN  LBS/SQ.INCH) 

2034  FORMAT ( 8X  *4H  LAM*6X*9H  1/4  LOAD,6X,9H  1 /^HLOA D » 6X , 9H  3/4  LOAD*6X*l 

10H  1  LOAD) 

2035  FORMAT ( 2  7X  * 26H  SHEAR  STRESS  TOTAL  IN  PSI) 

2036  FORMAT ( 9X *  I  3  *4F 1 5 .6) 

2037  FORMAT (9X,I3*4F15.2) 

2038  FORMAT ( 6X  * 14H  GAUGE  LINF  =  *I3*19H  LOAD  INCREMENT  =  *13) 

2039  FORMAT ( 33X  *  1 2H  (CORRECTED)) 

2040  FORMAT ( 9X  *  1 4H  COMPRESSION  =*E8.2*5H  KIPS*4X*10H  TENSION  =,F8.2*5H 
1 K I  PS  *4X  *  1 2H  SUM  HORIZ  =,F8.2,5H  KIPS*4X*12H  C/T  RATIO  =*F5.2} 

2  042  F 0 RMA T (8  X  ,4 H  L  A M  *  6 X  *  1 0 H  1  L 0  A  D  ) 

2044  FORMAT ( 9X* 13 ,F15.6 ) 

2045  FORMAT ( QX  ,  13 *F15.2  ) 

~TO$£~FORMAT(  12X,T0H~  INT.MOM  ^TfToT.2*5H  IN-K,3X,10H  EXT. MOM  =*F1o.2,5H  I 
1 N-K  *  1 2H  SUM  MOM  =»F10.2*5H  IN-K*4X*18H  IN/EX  MOM  RATIO  =*F5.2) 

2052  FORMAT ( 1X*51H  GAUGE  LOAD  COMPRESSION  TENSION  SUM  OF*6 

IX, 4H  C/T  *8X  *9H  INTERN AL  *6X  *9H  EXTERNAL *8X ,7H  SUM  OF,4X,8H  INT/EXTT 
2054  FORMAT ( 2X  *  1 2H  LINE  INCRE*5X,6H  FORCE, 6X,6H  FORCE, 9X,4H  T+C,6X*7H 
1  RATIO*8X,7H  MOMFNT*8X,7H  MOMENT, 8X,8H  MOMENTS, 5X ,6H  RATIO) 

2056  FORMAT ( 20X,4H  (K),8X,4H  (K),10X,4H  (K),21X,7H  (IN-K),8X,7H  (IN-K), 
18X,7H  (IN-K)) 

2  058  FORMAT  (  4X,I2*5X,U  *  1  X  ,  3F1 3 . 2 , 5  X  *  F5 . 2 , 2  X  ,  3F 1  5  •  2 , 5X  ,  E5  •  2  ) 

2  060”  FORMA  TT7X , 12 H  GA UGE  LINES ,13, 3H  TO,  I  3) 

2062  FORMAT ( 30X,20H  SHEAR  FORCE  IN  KIPS) 

2064  FORMAT ( 8 X , 35H  TOTAL  SHEAR  TOP/TOTAL  SHEAR  BOTTOM)  _ 

2066  FORMAT ( 19X,  I3,3F15.2) 

2067  FORMAT ( 8X  *  1 7H  SHEAR  GAUGE  LINE*6X,9H  1/4  LOAD,6X*9H  1/2  LOAD*6X,9H 
1  3/4  LOAD , 6 X , 1 OH  1  LOAD) 

2068  FORMAT ( 8X , 1 7H  SHEAR'GAUGE  LINE*6X, 1  OH  1  LUADT” 

2070  FORMAT ( 12X, 13, 3H  TO , I  3  ,4X * 4F 1 5 • 2 ) 

2072  FORMAT(  12X,I3,3H  TO  ,  I  3 ,4X  ,F  1 5  *2) _ 

2080  FORMAT(10X,50H  MODULUS  OF  ELASTICITY  VALUES  USED  IN  CALCULATIONS) 

2082  FORMAT ( 1  OX ,4H  LAM,7X,8H  F ( COMP ) , 7X * 8H  E ( TFNS )  ) 

2083  FORMAT  (  22X  ,6H  (PSI),9X,6H  (PSD)  _ _ 

2084  FORMAT  (  1 2X ,1 2 *6X ,d'0.0 ,5 X  , E10 •  0) 

2086  FORMAT (7X ,11H  GAUGE  LINE, 13, 3H  TO, 13) 

2090  FORMAT ( 35X ,60H  COMPARISON  OF  C/T  AND  INT/EXT  RATIOS  FOR  GAUGE  LINE 
IS  11  TO  ,13) 

2091  FORMAT ( 35X,61H  COMPARISON  OF  C/T  AND  INT/EXT  RATIOS  FOR  GAUGE  LINE 
IS  3  TO  16) 

2  092  FORMAT  (  2  5X“*  105H  LOAD  I NCRE •  LOAD  INCRE.  TO AD  I NCRE •  LOAD 

1  INCRE.  LOAD  INCRE.  LOAD  INCRE.  LOAD  INCRE.  ) 

2094  FORMAT (30X*2H  2*13X,2H  3*13X,2H  4*13X*2H  5*lOX,8H  2 +3 +4+5 ,8x,6H  3+ 
1 4  +  5 , 1  OX  * 4H  4+5) 

2096  FORMAT ( 2X  *  1 1 H  C/T  R A T I OS  * 6X * 7F 1 5 . 3 ) 

2098  FORMAT ( 2X, 15H  INT/EXT  R A T I  OS  ,  2X , 7F 1 5 • 3 ) 

3000  FORMA  )  (  10X, 39H  EQUATTPN5~PFFTN BFsr~ETT  OF  STRAINS* 

3002  FOR-MAT  (  1  5X  ,  1 7H  LOAD  INCREMENT  =  ,  1 3 , 5X  ,  E 1 2 .4 , 5H  +  F12#4,PH  Y*5H 

1  +  ,E12.4*5H  Y**2 )  _ 

3004  FORMAT ( 39X , 9H  (FITTED)) 

3006  FORMAT ( 6 IX , 9H  (FITTED)) 

4201  FORMAT ( 53X ,25H  1  INCH  =  600  LBS/SQ.INCH)  _ 

- ZTZU2 — h  URMA“T(  5X,1H.,9X*1H.,vx*1H.  ,9x*  ITT.  ,9XTTH.  VOX  ,1H.  ,9X  ,TH.T0')T5“TH.  *nX, 


AT-') 


X  >  .  1  A  I  \  -  ’  ,  X  (> 


^  . 

(  F  . 

30 

Xl\ 


•v\ 


tip 


I  j 


0.1 


..  ..  . ,  q  ,\  \  ,  3  |  Q  p  . 

I  x  «  TAG  . 

' i .  j  m  ao  h  r  r • x r ) t a m q (  - 

(  f  '*  1  *  I  ^  “  1  ’  91 A  1  *  c  ,  x  ‘  F  1  T  ‘  M  ci  p 

t  *  I  ^  I  A  vl  ^  Mrr,XPC)T4vq(  ' 
dn-'-,0'\^'  /  T  '  ?  3  9V '  P  MAC  »  X'  P  )  T  A  M  P  C 

HP  ♦  X ?  «  !  A O  i  •>•  \  i  <tV  ‘  '  J  :J  *  )  T  A  v 

f  0 AO  J  r  hoi 

(Tpq  [/||  JATQr  •'  •'  3  *  T  '•  flA^H?  H  ♦  X  T  S  )  T  A  Mfl  03 

^  r  *XP  )  T AM90 
(^.?f^tfI»VP)  TAMflO 

(  '.  HPX*fT*  «  3HIJ  30UA0  HA  f  #  Xd  )  T  A  M,90  3 

no3To?flfloot  h<^  r ♦  xc^ ) tah^oi 
•'•'•■  .  I)|  H?tS«  •  ’  '  -  HA r  •  x  -  )  TAMA03 

Hcr»n\,  'qj>  \i  -inH  mu  ’  h  s  r  •  x  a  *  -:*.  °  T  >  t 

( f!AOJ  T  H'  r  ♦  Xr) •  v  a  !  ha*  '  P  )  T A m p 0 1 

(  A,r  '  ''  )  TflVPAl 

(C.rrn*PT*X'-)  TAvqoO 
=  MO  •  T>  .■■•  ■.  .  .  .  .Tl  !  I  ftXSX)TAMfl03 

3 \ H I  H 8  X • X  A  «  X-H I  •  .  • *  r .  ■ - 

101  3£N  I  AO  ■  r  1 

...  ....  .  .  ’ 

.  ■  .  '  .  IdfXdt  *  '  r  •  >  )  A 

.  .  .  .  .  . 

K  ft  *  »V*XIS*  (X)  H  A  «  X  0 1  •  (  '"  '  HA*  X  • 

<(X-MI)  HTtXftI 

' .  ^  -  •  X  ?  •  S  •  £  I1C  ♦  XT#  •  •  X  11 

(fT,Cr  Hr.  c'  T  .  ^  ~A  T  J  "iP  '  '  Mil  *  Y' )  TAMP33 
(pqiX  / 1  POAOl  fl‘3u"  HOP  #  X0C  )  TAMP03 
'  '  •  - 

(  \  •  ?  r  ^  r  *  F I • X P I ) T  MPO 

•  OAOJ  A\f  H9*Xd*3HIJ  30UA0  flA3H  •  :  ' 

( gao j  r  Hor*Xc>*aAOj  a'f  r 

~'0  r  *  j  3-  !A0  W.  P  •  ■  HT  f  *  x  c  )  T  AM  V  1 1 

(  P.p  f  3,i.X»',FT*PT  HF#FT*XSf  1  T  A  <iqO"> 

(  p  .  ?  n  *  X  A  ,  C  T  .  C  r  ■  If  •  *  !  •  X S  X  )  T  Avqo  3 
T  '  JinjAO  M  T  C1?I!  ?3'  JAV  VTn!TAAJ“!  H)  5»lJt:00M  HOT*  X  P  [  )  T  A  M  q  O  ^ 

<  <  '  ‘  T  )  l  H°*  XT  .  (qMC  ))  HF*Xr*MAJ  H  A •  XO  f  )  T  AMH03 

Hd*XP*(T’M  HA*  XS5  >  TAW^C'l 
(  .  r  -j.  #  ;  “J*  XA*  :  *  <q  r  )  Tav-'O^ 

(fr,nT  Hr ,  c  j , -n*  t  j  -jp  ;a  o  h  r  r  *  xt  )  t  Avar  i 

I  201  T  A  '  *  \  T  1  y  OH  A  TV)  ^0  HO?.]  .  ■  ■ 

(f!#  o’  rr  ?i 

•  •  y T HI  0 H A  TND  30  HO? I  B  >  HI d* X?£ ) TAMA03 

(?f  OT  c  ?X 

.  .  .  )  f  *  X?  S I T  AMH03 


^ >  r 
^  H  ♦  = 

'At  |P 


<01 
[  )  Hr 


n  ao 


Y  *' 


.  -a  ',T 

. 


•  <  .Y  H< 


l  p  .  t  *  ; 

•  <:  n 


T  5 


^0 


•  #3f*^HI  OAOJ 

•  HS'Xei*A  H5#XfX*l  ^  r  •  V  J*  •  ,  s;  HS«X0£ITAM 

(  3  •*  HA  •  XP  r  ♦  r3  +  A  f 

.  1  v  •  v  -I .  P 01  T  Afl  T\D  HI  I« XS  ) TAM '-'r  1 
If«?I3T*XS*?0ITAfl  TX3VTHI  Hi  I#X 
r  j  -i  tp^o  OH  I  Ml  ’MO  I  T  AMi.v-i  hp«**X'  f)TA*  9:0^ 


4  •  9  A  •  * )  .  *  .  = 


nr  I  HI  r • X^  I  )  TAMS  ■ 
(q**y  +  X 

(  f  r">  T  id  HOfXf  "  )  T a v >! • ' ' ; 

(  (  03 T T  I  1 )  HP#  X  :  o  )  TAMflO^ 

HI  •  ) 2 \ ? fl  J  s  HOH  t  r  >“  J  '  . 


.  r 


•  • 


...  . . . . .  • 


•  • 


•  1  >  f  A  V  A  '  n 


dSO<; 
pf  n<; 

-  c 

^  p  n  r 
Ff  0^ 
AFOS 

?F0^ 

^FP^ 

VFDC 

rco<; 

npo<; 
0  A  P  ^ 

AAO^ 

PADC 

2APS 

A?P,^ 

8?r<; 

p^<; 
A?0  A 
VdO^ 

ftpPF 

OftOF; 

prp<; 

F8PC 
Aftp<; 
28  OF 
OP  PF 

IPOS 

SPOS 

APOS 

dPOF 

8P0S 

nOPF 

SOPf 

AOr  r 

IOSA 
crj  c  ■ 

J  j  r*  • 


C-15 


11H. *9X* 1H.  »9X*1H. *9X , 1H. >9X  » 1H. ) 

4203  FORMAT ( 5  0X  *31H  1  INCH  =  300  MICRO  INCHES/INCH) 

4204  FORMAT!  3X*127H . . 

]• • . . .  . ... . . .... . .  .  .  ...  .  . . .  .  . 

2.  . . ) 

4206  FORMAT ( 50X 1 18H  1  INCH  =  100  P SI») 

4208  FORMAT! 11X*17H  LOAD  I NCRFMENT  * ' *  12) 

4246  FORMAT  ! 4SX  *  1 H •  ) 

4266  FORMAT ( 63X  *  1 H. ) 

°  00 1  FORMAT ( 1 X  * 1H  ) 

9002  FORMAT(IHJ) 

9003  FORMAT !  1HK  ) 

9010  FORMAT ( 1 H 1 ) 

STOP 

END 
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PART  II  SUBROUTINES 


SUBROUTINES  USED  ARE  AS  FOLLOWS. 

1.  STRAIN 

2.  STRESS 


3.  FORCE 

4.  MOMENT 

5.  COMPAR 

6.  FIT 

7.  PLOT 

8.  DELTA 

9.  CORREC 
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$  I BFTC  NIARTS 

SUBROUTINE  STRAIN 
C 

C  T H IS  SUB ROUT  I  NEC A LCUL A TES  T H F  S T R A  I N  0 F  EACH  GAUGE  L E NGTH  FROM 

C  THE  TEST  DATA#  SUBROUTINE  STRAIN  ALWAYS  RETURNS  THE  PROGRAM  TO 

C  THE  TEST  STRAINS  AND  IS  INDEPENDENT  OF  ALL  OTHER  SUBROUTINES* 

C 

C  REQUIRED  IN  COMMON  DR( 30*20*5) *ST (30*20*5) 

C 

COMMON  DR ( 30  *  20  *  5 )  *ST ( 30  * 20*5) *DEL ( 24*20*5 )  *DEL T (24  *  20*5 T * 

1 HORC (30*5) *HORT(^0*5)  *SHOR(30,5)  ,RHOR(30*5)  * TE I MO ( 30  * 5 ) , F I MO ( ? 0 ) * 
1 FXMO  (30*5)  *  SMOM  ( 30  *  5  )  ,RM0M(30*5)  *EC(20)*FT(2^>)*P(?0)*A(20)*Y(20)* 
lX(3r')  ,  BB 1  (30*5)  *BB2(30*5)  ,BB3  (30*5)*C1(30*5)*C?(30*^),C?(^0*5)* 
1AA1 >AA2* AA3  *AA4*AA5  *AP2*AP3  >AP4  *AP5  *BP1 *BP2  >BP3  * 

1SH2  *SH3*SH4*SH5  *SM2*SM3*SM4,SM5  *TH1 ,TH2 , TH3  * TH4 , TH5  * TH6 ,TH7 , TM 1 , 

1 TM2  *  TM3  * TM4 *  TM5 *  TM6  *TM7  *DI V  *  B1 * B2  *B3  *DET * 

1B*D*CAL  *TL*GAGL ♦ LIMIT  *LOLIM*LAML*LSP AN* LOAD *PLOAD* SPAN* NGL  *NLAM, 
1NL,NNL*NGLM*NGLG,NGLL*NGLR,R,I *  I N ,NNN , NORUN 


DO  310  NL  =  2  *NNL 
DO  312  NGL  =  1  * NGLM 
DO  314  NLAM=1 *20 
ST(NGL*NLAM,NL)=0.0 
314  CONTINUE 
312  CONTINUE 
310  CONTINUE 
NLL  =  1 

DO  308  NL  =  2  *NNL 
DO  306  NGL= 1 *NGLM 
DO  304  NLAM=1 ,20 

ST  (NGL,NLAM,NL  )  =  (DR (NGL*NLAM,NL)-DR (NGL*NLAM,NLL )  )*CAL 
304  CONTINUE 
306  CONTINUE 
308  CONTINUE 
RETURN 
END 


S I BFTC  SSERTS 

SUBROUTINE  STRESS 

C _ _ _ 

C  THIS  SUBROUTINE  CALCULATES  THE  STRESS  OF  EACH  GAUGE  LENGTH. 

C  SUBROUTINE  STRAIN  OR  SUBROUTINE  FIT  PRECEEDS  SUBROUTINE  STRESS. 

C  _ _ _ _ _ _ 

Z~  REQUIRED  in  COMMON  DR(30*20*5  ) *ST(30,20*5) *EC(20)  ,ET(20)  *P(20) 

c 

COMMON  DR(30*20*5) *ST(30*20*5) *DFL(24*20*5) ,DELT(?4*70»5) * 

1H0RC (30*5) *HORT ( 30*5 )  ,SHOR ( 30,5 )  *Rh0R( 30*5 ) * TF I MO ( 30  *  5 ) ,FIMO(20) * 
1EXM0( 30*5 ) *SMOM( 30*5 )  *RMOM( 30*5 )  ,EC ( 20  > *ET ( 20 ) ,P ( 20 ) * A ( 20 ) *Y ( 20 ) * 

1X130) *BB1 (30*5 ) *RR2 ( 3  0*5) *BB3(3  0*5 )  ,C1 (30*5 ) *C2( 30*5 ) , C  3 ( ^0*5  )j _ 

1AA1 *AA2*AA3*AA4*AA5*AP2*AP3*AP4*AP5*BP1*BP2  ,BP3* 

1SH2,SH3*  SH4  *  SH5 , SM2 , SM3  *  SM4 , SM5 , TH 1 , TH2 , TH3 * TH4 , TH5 , TH6 , TH7 , TM 1 , 
1TM2  *TM3 *TM4*TM5 * TM6  * TM7 * D I  V  * B1 *B2>B3*DET  * 

- IB  *D»CAL  *  TL  * G7TG’ C"*TTI^TT~ ►TTrrTTrnrA'MrTET^'ffTNTirOTCBT^^ 
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D 

.  ('■'^c,  sc  )T  J-30,  (  p .  r  c  .  a  S  )  IT*  (  ?*  cp,  Or)  jp  ,  (  P*OS»OF  )^fl  IAOMMO*) 

'  M 1 3 # (  # OF )  M I IT # ( P# OF ) AOhA#  t  A* OF ) A0H2 1 ( ?• 0^ ) T AOH# (^ • Of ) 3A0HI 
•  (  0^  )  A#  (  os  )  T3#  <  0?  )  33#  c  .  ■  c  #  (  #0f  )  vn  <!£#<?# OF  >0MX3f 

.  (  )c*0,  (  %OF  'CO.  (  ^,0f)  ro.  (  ?  ♦  O  C  )  ►'  -!0  ,  (  p  #  0  ■  )  ^  PP  *  (F.OF)inP#  (  o  F  )  X  I 

.  q  . c  #  far*  *  P  r  #  a  o  A  #  c  n ' • S  q  A  #  P  A  A • A  A  A • F  A  A • S  A  A •  fAAI 

IT#  ’•  IT  •  FHT  *  SHT  *  T  H  T  #  M  ?  ♦  A  M  '  •  ....  .  •  £ 

.  Tin.  c  *  r  r  »  V  I  ru  VMT  •  V/T  •  .?MT  •  AMU  C^T  »  ^mt  r 
« .)  .#  ■  -  - .  a  *  'iAo.it  ’a^  j...  (jj#  ’  i  r  » jo'r#  jr#  't  'f  .  i 


c 


1NL,NNL*NGLM,NGLG*NGLL »NGLR  »R> I  *  I N ,NNN , NORUN 


NLL  =  1 

DO  362  NL -2 *NNL 

DO  360  NGL  =  1  * NGLM 
DO  358  NLAM= 1 *20 
P(NLAM)=ST(NGL,NLAM,NL) 

IF (P (NLAM) )  352  *352  *  354 
352  ST(NGL*NLAM,NL)=P(NLAM)*EC(NLAM) 
GO  TO  356 

354  ST (NGL,NLAM,NL )=P( NLAM) *ET (NLAM) 
356  CONTINUE 
358  CONTINUE 


360  CONTINUE 
362  CONTINUE 
RETURN 

END 


$  I BFTC  ECROF 

SUBROUTINE  FORCE 
C 

C  THIS  SUBROUTINE  CALCULATES  THE  HORIZONTAL  FORCE  OF  EACH  LAMINATION 

“C  AT  EACH  GAUGEHLINET-  SUBROUTINE  STRESS”PRECEEDS  SUBROUTINE  FORCE. 

C 

C  REQUIRED  IN  COMMON  PR ( 30 »20»5 >  tST ( 30 *20 »5 ) *H0RC<30«5>  » 

C  HORT ( 30*5) *SH0R(30*5 ) , RHOR ( 30 , 5 ) * P ( 20 > * 

C  A (20) *EC(20) ,ET(20> 

C 

COMMON  DR (30  * 2  0  *5) *ST (30*20*5 ) *  DEL (24  *  2  0*5 ) *DELT (24*2  0 ♦ 5  >* 

1H0RC (30*5) * HORT (30*5 )  *SHOR( 30*5 ) *RHOR( 30*5 ) * TFI MO ( 30 *5  > *F I  MO ( 20 ) * 

1 EXMO (30*5) ♦ SMOM (30*5 > ,RM0M(30,5) *EC ( 20 ) ,ET ( 20 ) ,P ( 20 ) *  A ( 20 ) *Y ( 20 ) , 
IX (30) *BB1 (30*5) *BB2 (30*5 ) *BB3 (30,5 ) ,C 1 ( 30 *5 ) * C2 ( 30 ♦ 5 ) *C3 ( 30 ♦ 5 ) * 

1 AA1 * AA2  *  AA3 *  AA4* AA5  *AP2 * AP3 * AP4 * AP5  *RP 1 *BP2  *BP3* 

1 SH2  *  SH3  *  SH4 ,SH5*SM2*SM3*  SM4  *  SM5  *  TH1 * TH2 * TH3 , TH4 * TH5 * TH6 , TH7 , TM 1 , 
1TM2*TM3*TM4,TM5*TM6  *  TM7  * DI V  *B1*B2*B3*DET ♦ 

1B*D*CAL*TL*GAGL*LIMIT*L0LIM*LAML*LSPAN*  LOAD  *P LOAD  * SPAN  *NGL*NLAM* 
1NL*NNL*NGLM,NGLG*NGLL*NGLR*R»I*IN  *NNN  *  NORUN  __________  ________ 

C 

DO  480  NLAM= 1  *  20 

IF (NLAM.EQ.10.OR.NLAM.EQ.il )G0  TO  482 
A (NLAM ) =B*TL 

GO  TO  484 

482  A ( NLAM) =B*TL/2.0  _ 

484  CONTINUE 
480  CONTINUE 

DO  487  NL  =  2  *NNL 
DO  489  NGL= 1 *NGLM 

DO  485  NL AM= 1  * 20 

HORC (NGL*NL )*0.0  _ 

HORT (NGL*NL)=0.0 
SHOR (NGL*NL)=0.0 
RHOR (NGL  *NL )=0.0 
- 48b  CONIINUE - 
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489  CONTINUE 

487  CONTINUE 
NLL  =  1 

DO  486  NL  =  2  *NNL 
DO  488  NGL=1*NGLM 
DO  490  NLAM= 1  *  20 

ST (NGL,NLAM,NL )=ST (NGL ,NLAM,NL)*A(NLAM ) / 1000.0 

490  CONTINUE 

488  CONTINUE 
486  CONTINUE 

DO  470  NL  =  2  *NNL 

DO  472  NGL=1*NGLM  _ _  __  _ 

DO  474  NL AM* 1  *  20 
IF (ST (NGL*NLAM*NL) >476*476,478 
476  HORC ( NGL ,NL > =HORC ( NGL *NL ) +ST ( NGL ,NLAM , NL ) 

GO  TO  479 

478  HORT(NGL*NL)=HORT(NGL  *NL ) +ST ( NGL ,NLAM*NL ) 

479  CONTINUE  _ _  _ 

474  CONTINUE 

472  CONTINUE 
4.70  CONTINUE 

DO  481  NL=2*NNL 
DO  483  NGL  =  1 *NGLM 

SHOR ( NGL  *NL ) =HORT ( NGL *NL ) +HORC ( NGL , NL ) 

RHOR ( NGL *NL)=ABS( HORC (NGL,NL) /HORT ( NGL *NL ) ) 

483  CONTINUE 
481  CONTINUE 
RETURN 
END 


SIBFTC  TNEMOM 

SUBROUTINE  MOMENT 
C 

C  THIS  SUBROUTINE  CHECKS  THE  INTERNAL-EXTERNAL  MOMENT  CONDITIONS. 

C  SUBROUTINE  FORCE  PRECEEDS  SUB ROUT I NE  MOMENT. 

C 

c  REQUIRED  IN  COMMON  ST ( 30 , 20 *5  ) * Y ( 2 0 )  ,F IMO ( 20 )  , TF I  MO ( 3 0  *  5 > * 

C  X ( 3 0 ) *  EXMO (30*5 >  ,SMOM(30*5  > ,RMOM(30*5  > 

C 

COMMON  DR (30*20*5 ) *ST (30*20»5 ) *DEL ( 24  *  20 *5 ) »DELT ( 24*20*5 > * 

1H0RC (30*5) *HORT(30*5)  *SHOR ( 30*5  >  *RHOR( 30*5  > ,TFIMO( 30*5 > ,FIM0(20> * 
1EXM0 (30*5) *  SMOM (30*5 )  ,RMOM( 30*5 )  *EC ( 20  > * ET ( 20 ) , P ( 20 ) *  A ( 20 >  *Y ( 20 > , 
IX ( 30 ) *BB1 ( 30  *5 ) *RB2(30*5> *BB3(30*5 ) ,C1(30*5>*C2(30*5)*C3(^0*5), 
1AA1*AA2*AA3*AA4*AA5*AP2*AP3*AP4*AP5*RP1*BP2*BP3* 

1 SH2 , SH3  *  SH4 ,SH5*SM2*SM3,SM4*  SM5 , TH1 * TH2 * TH3 * TH4 , TH5 , TH6 , TH7 , TM 1 , 

1TM2*TM3,  TM4 , TM5  *  TM6  ,TM7  *DIV,B1 *B2*B3*  DET  * _ 

“TB»D»CAL*TL*GAGL*LIMITTCOLIM,LaML  *  L  SP  AN7LUAD  *P  LOAD  *SPANTNGL7N  LAM, 
1NL*NNL*NGLM*NGLG,NGLL  *NGLR*R  ,  I  ,  I  N  *  NNN  *  NORljN 
C 

NL AM=1 

Y ( NLAM )=D/2.0-TL/2.0 
DO  412  NLAM=2 ♦ 9 

—  YTNt  AMT* Y(NLAM-1 1-TL  _ 
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412  CONTINUE 
NLAM=10 

Y ( NLAM) =Y (NLAM-1 )-0 .75*TL 
NLAM=ll 

Y(NLAM) =TL/4.0 
NL AM=1  2 

~Y(NLAM)=Y(NLAM-1)+0.75*TL 
DO  414  NLAM= 1 3  »  2 0 
Y ( NLAM ) =Y( NLAM-1 >+TL 
414  CONTINUE 

DO  423  NL  =  2  *NNL 
DO  424  NGL  =  1  *  NGLL  *  2 
FIMO(NLAM)=0.0 
TEIMO(NGL*NL) =0.0 
DO  416  NLAM  =  1 *  10 

“FT  MO  ( NLAM ) =+ST ( NGL  ♦NO^RUT^YTML  AW  )'" 
416  CONTINUE 

DO  418  NLAM= 1 1*20 

FTBO'THL AW )  T T  N GL  *  N L  A M  ♦  N L ) * Y  ( N L  AM  ) 

418  CONTINUE 

DO  425  NLAM  =  1  *  20 

“TFTMO  ( NGL  *NL) *TF I MO ( NGL » NL ) +F I  MOJNUAM } 

425  CONTINUE 
424  CONTINUE 

DO  426  NGL  =  2  *  NGLR  * 7 
F I  MO (NLAM)=0.0 
TFIMO(NGL*NL>=0.0 
DO  419  NLAM=1*10 

FIMO (NLAM) =-ST (NGL  *NL AM ,NL ) *Y ( NLAM ) 

419  CONTINUE  _ 

DO  420  NLAM  = 1 1 *20 

FIMO ( NLAM ) =+ST ( NGL ♦ NLAM* NL ) *Y(NLAM ) 

420  CONTINUE 

DO  427  NLAM =1,20 

TFIMO(NGL*NL)*TFIMO(NGL*NL)+FIMO(NLAM) 
427  CONTINUE 

426  CONTINUE 
423  CONTINUE 

NGL  =  1 

X ( NGL ) =GAGL 

DO  430  NGL=3*NGLL*2 

X ( NGL ) =X ( NGL-2 ) +GAGL  _ 

430  CONTINUE 
NGL  =  2 

X ( NGL ) =GAGL 
“DO  432  NGL=4 ♦ NGLR  *  2 

X ( NGL  )  =  X ( NGL-2 ) +GAGL 
432  CONTINUE 

IE ( LOAD.EQ. 1 ) GO  TO  434 
IF(L0AD.EQ.2)G0  TO  436 
434  CONTINUE 

DO  438  NL=2.NNL 

IF (NL.EO.2)R=0.25*PLOAD 

IF(NL.EQ.3)R=0.50*PL0AD  _ 

IF (NL.EQ.4)R=0.75*PLOAD 
IE (NL.EQ.5 )R=1 .00*PLOAD 
DO  440  NGL  =  1  *  NGLL  *  2 
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440  CONTINUE 

DO  442  NGL=2*NGLR*2 

EX  MO  (NGL*NL):=‘-R/2*0*  (SPAN/2.  O-X(NGL) ) 

442  CONTINUE 
438  CONTINUE 
GO  TO  433 
436  CONTINUE 

DO  444  NL  =  2  *NNL 

IF ( NL*EQ*2 ) R=0.25*PLOAD 

IF (NL.EQ.3 ) R  =  0*  50*PL0 AD 

IF (NL.EQ.4)R=0.75*PL0AD 
IF (NL.EQ.5 )R=1 .00*PL0AD 
DO  446  NGL  =  1 *NGLL  *  2 

EXMO(NGL,NL)=R*(SPAN/?#0-X(NGL) )-R*(0.25*SPAN-X(NGL) ) 
446  CONTINUE 

DO  448  NGL  =  2  *  NGLR  *  2 

EXMO(NGL,NL)=-R*(SPAN/2*0-X (NGL) )+R*(0.25*SPAN-X (NGL) ) 
448  CONTINUE 
444  CONTINUE 
433  CONTINUE 

DO  450  NL=2*NNL 
DO  452  NGL  =  1  * NGLG 

SMOM(NGL*NL)=TFIMO(NGL*NL)+EXMO(NGL,NL ) 

RMOM (NGL*NL)=ABS(TFI MO ( NGL  *NL ) /EXMO ( NGL  *  NL ) ) 

452  CONTINUE 
450  CONTINUE 
422  CONTINUE 
RETURN 
END 
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SUBROUTINE  COMPAR 

THIS  SUBROUTINE  SUMS  THE  C/T  AND  INT/EXT  RATIOS.  THE  GAUGE  LINES 
C  SUMMED  WERE  DETERMINED  FROM  A  NON  INTERFERENCE  OF  LOAD 

C  CONCENTRATION  AT  THE  WORKING  STRENGTH  LOAD  LEVEL.  THEREFORE  GAUGE 

C  LINES  11  TO  16*  20*  OR  24  WERE  USED  FOR  THE  SINGLE  LOAD  CASE  AND 

C  3  TO  16  WERE  USED  FOR  THE  TWO  POINT  LOAD  CASE.  SUBROUTINES  FORCE 

C  MOMENT  MUST  HAVE  BEEN  CALCULATED  PREVIOUSLY  IN  THE  PROGRAM. 

C _ _ _ 

C  REQUIRED  IN  COMMON  RHOR ( 30*5 ) *RMOM( 30*5 > 

C 

COMMON  DR (30*20 *5)  *ST(30*20*5) » DEL ( 24  *  20 *5 ) *DELT ( 24 ♦ 20  * 5 >  * 

1 HORC (30*5) *HORT(30*5>  *SHOR(30*5) ,RhOR(30*5)  *TFIMO(30*5>*FIMO(20)* 
1 FXMO (30*5) *SMOM(30*5>  .RMOM (30*5) *  EC ( 20 > ♦ ET ( 20 > . P ( 20 ) *  A ( 20 ) . Y ( 20 > * 
1 X ( 3 0 ) *BB1 (30*5 ) *BB2 (30*5 ) *BB3 (30*5 ) »C1 ( 30*5  > »C2 ( 30*5  > »C3( 30*5 )  » 

1 A a 1  *  AA2 *AA3*AA4*AA5*AP2*AP3*AP4*AP5  *BP1 *BP2  *BP3  * 

1SH2  *SH3  *SH4  *SH5  *  SM2  *SM3  *SM4  *  SM5  *  TH1 ,TH2 , TH3 , TH4 , TH5 * TH6 , TH7 , TM 1  * 
lTM2»TM3*TM4*TM5*TM6»TM7*DIV»Bl*B2»B3»DETt 

1B*D*CAL*TL  *GAGL  *LI M I T  *  LOL I M  *  LaML  *  L SPAN  *  LOAD  *P LOAD  *SP AN *NGL*NLaM* 
1NL*NNL*NGLM*NGLG*NGLL,NGLR,R*I ,IN*NNN,NORuN 
C 
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SM?=0.0 
SH3=0*0 
SM3=0 • 0 
SH4  =  0  #  0 

SM4=0* 0 
SH5=0*0 
SM5=0.0 
TH!=0*0 
TM 1  =0  • 0 
TH?=n.O 
TM?=0 • 0 
TH3  =  0  *0 
TM3=0*0 
TH4=0*0 
TM4=0.0 
TH5=0.0 
TM5=0.0 
TH6=0.0 
TM6=0.0 

TH7=0*0 
TM7  =  0«  0 

“TFT  LOAD • F Q • T ) GO  TO  518 
IF (L0AD*EQ«2 ) GO  TO  519 

518  CONTINUE 
NL  =  2 

DO  520  NGL=11*NGLG 
5H2=RHOR(NGL*NL )+SH2 
SM  2  =  RMOM ( NGL  *  NL )+SM2 
520  CONTINUE 
NL  =  3 

DO  522  NGL=11*NGLG 
SH3=RHOR(NGL»NL)+SH3 
5M3=RMOM(NGL*NL)+SM3 
522  CONTINUE 
NL-4 

DO  524  NGL=11*NGLG 
c,H4  =  RH0R  ( NGL  *NL  )  +SH4 
SM4=RM0M ( NGL ♦ NL ) +5M4 
524  CONTINUE 
NL  =  5 

DO  526  NGL=1 1 *NGLG 
SH5=RH0R (NGL»NL )+SH5 
SM5=RM0M(NGL*NL )+SM5 
526  CONTINUE 

IF(NGLG.EQ.24)DIV=14#0 

IF(NGLG.EQ*16)DIV«6.0 

IF(NGLG.E0.20)DIV*10.0 
GO  TO  538 

519  CONTINUE 
NL  =  2 

DO  530  NGL  =  3  *16 
SH2*RH0R (NGL *NL)+SH2 
SM2=RM0M ( NGL  *NL ) +SM2 
530  CONTINUE 
NL  =  3 

DO  532  NGL  =  3  *  1 6 
SH3=RH0R(NGL*NL )+SH3 
- 5  M  3  = RMUTTT7TG  LT)TnT?OTg~ 
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532  CONTINUE 
NL  =4 

DO  534  NGL  =  3  *16 
SH4  =  RH0R  T NGL , NL )+SH4 

SM4=RM0M ( NGL ,NL )+$M4 
534  CONTINUE 
NL  =  5 

DO  536  NGL  =  3  *16 
SH5=RH0R (NGL ,NL )+SH5 
SM5=RM0M ( NGL  *  NL ) +$M5 

536  CONTINUE 
D I V=14.0 
538  CONTINUE 

TH1=SH2/DIV 
TM1=SM2/DIV 
TH2=SH3/D IV 

TM?=SM3/DI V 
TH3=SH4/DIV 
TM3=SM4/DIV 
TH4=SH5/DIV 
TM4=SM5/DI V 

TH5= ( SH 2+SH3+SH4+SH5 ) / (DIV*4.0) 
TM5  = ( SM2+SM3+SM4+SM5 ) /(DIV*4.0) 
TH6=(SH3+SH4+SH5 ) / ( DI V*3.0> 

TM6= ( SM3+SM4+SM5 ) / ( D I V*3 • 0 ) 

TH7  =  ( SH4+SH5 ) /(DIV*2.0) 

TM7= ( SM4+SM5 )/(DIV*2«0) 

RETURN 

END 


SIBFTC  TIF 

SUBROUTINE  FIT 


C 

c 

c 

c 

c 

c 

c 

c 

c 


c 


THIS  SUBROUTINE  CALCULATES  THE  BEST  FIT  SECOND  DEGREE  CURVE 
THROUGH  THE  STRAIN  DATA.  THEN  THE  STRAIN  VALUES  ARE  REVALUATED 

AS  A  FUNCTION  OF  THE  BEST  FIT  CURVE.  SUBROUTINE  STRESS  PRECEEDS 
SUBROUTINE  FIT. 


REQUIRED  IN  COMMON  Y ( 20 ) * BB 1 ( 30 , 5 )  * BB2 ( 30 ,5 )  ,BB3 ( 30 » 5 ) *P ( 20 ) , 

ST ( 30 ,20,5  > ,C1 ( 30,5 ) ,C2( 30,5 ) ,C3(3o,5> _ 


COMMON  DR (30*20*5), ST (30, 20*5) , DEL (24*20,5 )  ,DELT ( 24 , 20 ,5  >  , 

1H0RC (30,5) ,HORT (30, 5 )  ,SHOR ( 30,5 ) ,RhOR(30,5 ) , TF I  MO ( 3 0 , 5 ) , F I  MO ( 2 0 ) * 
1 FXMO (30,5) ,SM0MC*0,5> ,RM0M(30,5) ,EC ( 20 > , FT ( 20 > ,P ( 20 ) , A ( 20 ) , Y ( 20 ) * 
1X(30) »BB 1 (30*5) *RB2(30,5) *BB3 (30*5 ),C1(30*5)*C2(30*5>*C3(^0,5), 

1 AA1 *AA2  * AA3  *AA4*AA5  *AP2  *AP3  * AP4  * AP5  *BP 1 *BP2  *BP3  * 

1 SH2  * SH3  * SH4 , SH5 ♦ SM2  * SM3  *SM4  ,SM5  *  TH1 ,TH2  *TH3 , TH4 , TH5 , TH6 , TH7 , TM1 , 

1TM2*TM3*TM4*TM5*TM6,TM7*DIV*B1 *B2*B3*DFT , 

1B*D*CAL  *TL*GAGL* LIMIT  *L0LIM»L AML ,LSPAN , LOAD *PL0AD* SPAN, NGL ,NLAM, 

1NL,NNL,NGLM,NGLG,NGLL ,NGLR,R , I , IN, NNN, NORuN 
DO  800  NL AM=1 ,20 

Y/NLAM I  =(->.■(■) - — —  . . . . 
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800  CONTINUE 
NL AM=1 

Y(NLAM) =D-(TL/2«0) 

DO  802  NL AM =?, 9 
Y(NLAM) =Y(NLAM-1 )-TL 
802  CONTINUE 
NL AM=1 0 

Y(NLAM)=Y(NLAM-1 )-0.75*TL 
NLAM^ll 

Y ( NL AM ) =Y ( NL AM- 1 ) -0  *8  0*f  L 

NLAM=12 

Y(NLAM) =Y(NLAM-3)-2.0*TL 
DO  804  NLAM= 1 3*20 
Y(NLAM) =Y(NLAM-1 )-TL 
804  CONTINUE 
AA1 =0  •  0 

AA2=0.0 
A A3=0 • 0 
AA4  =  0  •  0 

AA5=0.0 

AA1=20.0 

DO  806  NLAM=1 9  20 

AP  2  =  Y ( NLAM ) 

AP3=Y(NLAM)**2 
AP4=Y(NLAM)**3 
AP  5  =  Y (NLAM ) **4 
AA2=AA2+AP2 
AA3=AA3+AP3 

AA4= AA4+ AP4 
A A5=AA5+AP5 
806  CONTINUE 

DET=AA1 * ( ( AA3*AA3)~( AA4*AA4 ) )-AA2*< ( A A2*AA5 ) - ( A A3* A A4 ) )+AA3*( (  AA2* 
1AA4)-(AA3*AA3 ) ) 

DO  808  NL=1 *NNL 

DO  810  NGL  =  1 *NGLM 

BB1 (NGL*NL)=0.0  _ __ _ _  _ 

BB2  (NGL*NL)=0«0 
BB3(NGL*NL)=0.0 
810  CONTINUE 
808  CONTINUE 

DO  812  NL  =  2  *NNL 
DO  814  NGL=1*NGLM 
DO  816  NLAM  =  1 *20 
P ( NL AM ) =ST ( NGL *NL AM  *NL ) 

BP  1 =P ( NL AM ) 

BP  2=P ( NLAM ) *Y (NLAM ) 

BP3=P(NLAM)*Y(NLAM)**2 
BB1 ( NGL  *NL ) =BB1 (NGL*NL)+BP1 
BB2 (NGL*NL)=BB2(NGL*NL)+RP2 
BB3 (NGL  *NL ) =BB3 ( NGL  *NL ) +BP3 
816  CONTINUE 
814  CONTINUE 

812  CONTINUE 

DO  818  NL  =  2  *NNL 
DO  820  NGL=1*NGLM 
BlaBBl ( NGLtNL ) 

B2  =  BB2 ( NGL  *NL ) 

- __ — y  3  ^  ^  CNOL  VNTX .  . . . . . 
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^#(mAJM)Y*(VA  IM)  q=rqn 

r r,n+(  jm.  jr  m  rsp=  (  jm.  jom)  r 

iM.  0M)^aqa( jm* joM)^qq 
(  jm. jom ) <  jm*  jom) pqp 

qijMiTMoo  ^tr 
q  MTTl/OO  A  r  ° 
qilMITMOO  ^r« 
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VJOM# r= JOM  n^R  on 
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C-25 

Cl  ( NGL  *  NL )  =  ( R 1* ( (AA3*AA5)-( AA4*AA4) )-AA2*(  ( B2*A A5 ) - ( B3*A A4 )  )+AA3* 
1 ( (B2*AA4)-(B3*AA3) ) )/DFT 

C2(NGL*NL)=( AA1*(  ( B7*  AA5  ) - ( B3*A A4 )  )  -B 1  *  (  (  A  A  2*  AA5  )  -  (  AA7*A  A4  )  )+AA7* 
1 ( (A  A  2 *  B  3  )~( A  A  3 *  B 3 )  ) ) /OF T 

C3 (NGL*NL)*( AA1*( ( AA3*B3 ) -( AA4*B2 ) )-AA2*( (AA2*B3 )-( AA5*B2) )+Bl* 

1 ( ( AA2*AA4)-( AA3*AA3 ) > ) /DFT 
820  CONTINUE 
818  CONTINUE 

DO  822  NL  =  2  *NNL 
DO  824  NGL=1 tNGLM 
DO  826  NLAM  =  1 *  20 

ST (NGL,NLAM,NL )=C1 ( NGL ,NL ) +C2 ( NGL ,NL ) *Y ( NLAM ) +C3 ( NGL ,NL ) * ( Y ( NLAM ) 
1**2) 

826  CONTINUE 
824  CONTINUE 
822  CONTINUE 
RETURN 
END 


SIBETC  TOLP 

SUBROUTINE  PLOT 


C _ _ 

C 

C  THIS  SUBROUTINE  PLOTS  STRAIN,  STRESS  OR  SHEAR  STRESS  DISTRIBUTION. 

C  THE  VALUES  ARE  DETERMINED  IN  THE  MAIN  PROGRAM.  SUBROUTINE  STRAIN 

C  OR  SUBROUTINE  STRESS  OR  SUBROUTINE  DELTA  MUST  PRECEED  SUBROUTINE 

C  PLOT. 

C  REQUIRED  IN  COMMON  P(20) 

“C~ 

COMMON  DR  (30  *2-0  *5)  *ST  (30*20*5  )  *DEL(24*20*5  )  *DELT  (  ?4  *  20  *5  >  * 

1HORC (30*3 ) ,HORT( 30*5 ) »SHOR ( 30*5 ) ,RHOR (30,5)*TFIM0(30*5>,FIM0(20), 

1 EXMO (30*5) *SM0M(30*5)  ,RMOM(?0*5)  *EC  (  20  )  * ET ( ?0 ) * P ( 20 ) , A ( 20 )  , Y ( 20 > * 
IX ( 30 )*BB 1(30*5)* BB 2(30*5) *BB3 (30,5 )  ,C 1 ( 30 , 5  > , C2 ( 30 , 5  > ,C3 ( 30 , 5 )  , 

1 A A 1 ♦ AA2 *AA3  *AA4*AA5  *AP2*AP3  »AP4  *AP5 ,BP1 *BP2  *BP3  * 

1 SH2 , SH3 , SH4 , SH5 ,SM2,SM3 ,SM4,SM5 , TH1, TH2 , TH3 , TH4 , TH5 , TH6 , TH7 ,TM T , 
1TM2*TM3  *  TM4 ,TM5,TM6,TM7*DIV,B1 *B2  *B3,DET , 

1B*D*CAL ,TL ,GAGL *  LIMIT  *L0LIM,LAML ,LSPAN ,L0AD*PL0AD* SPAN ,NGL  *NLAM, 
1NL  ,NNL,NGLM,NGLG,NGLL  ,NGLR,R  ,  I  ,  I  N  ,  NNN  ,  NORijN 


IE (P (NLAM) .LE.1815 .0. AND. P( NLAM) .GT.1785 
IE (P( NLAM) .LE.1T85.0. AND.P( NLAM) .GT.1755 
IF (P (NLAM) . LE. 1 755. 0. AND. P( NLAM) .GT.1725 
IF(P(NLAM) .LE. 1725. 0. AND. P( NLAM) .GT.1695 
IF (P( NLAM) .LE. 1695.0. AND.P( NLAM) .GT.1665 
IF (P (NLAM) .LE. 1665.0. AND.Pt  NLAM) .GT.1635 
IF (P (NLAM) .LE. 16 35.0. AND.P( NLAM) .GT. 1605 
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IF(P(NLAM)  .LE.  139*5.0.  AND.P(NLAM)  .GT.  1365. 0)WR  I TE  (  6  *4020) 

I F ( P ( NLAM ) .LE.1365.0.AND.P( NLAM) .GT . 1 ^35 .0 ) WR T TE ( 6  *  40? 1  ) 

IF (P (NLAM) .LE. 1335.0. AND.P ( NLAM) .GT.l 305.0) WRITE (6 *402 2 ) 

IF ( P ( NLAM) .LE. 1305 .0. AND. P (NLAM ).GT . 1 275 .0 ) WR I TE (6 *4023 ) 

IF(P(NLAM) .LE. 1275.0. AND. P( NLAM) .GT. 1245.0) WRITE (6 *4024 ) 

IF (P (NLAM) .LE. 1245 .O.AND.P ( NLAM) .GT . 1 ? 1 5 .0 >WR ITE ( 6 *4025 > 

IF (P (NLAM) .LE. 12 15.0. AND. P< NLAM) .GT. 11 35.0 )WRTTF (6*4026 ) 

IF (P( NLAM) .LE. 11 85 .O.AND.P (NLAM) .GT . 1 1 55 .0 ) WR ITE ( 6  *  4027 ) 

IF (P( NLAM) .LE. 11 56.0. AND. P( NLAM) .GT • 1 1 2 5 .0 ) WR I TE ( 6 ♦ 402 8 ) 

IE ( P ( NLAM ) .LE. 1 1 25 .0. AND. P ( NLAM) .GT. 1 095.0) WR I TE ( 6  * 4029 ) 

IF (P( NLAM) .LE. 1095 .O.AND.P (NLAM) .GT . 1065 .0 ) WR ITE ( 6 ♦ 4030 ) 

IF (P (NLAM) .LE.1065 .0. AND.P(NLAM) . G T. 1035. 0)WRITE(6* 4031) 

IF (P (NLAM) .LE. 1035 .0. AND.P ( NLAM).GT.1005.0)WRITF(6*4032) 

IF (P (NLAM) .LE. 1005 .0. AND.P ( NL AM ) . GT .0975 .0 ) WR I TF ( 6  *  4033 ) 

IF (P (NLAM) .LE. 0975 .O.AND.P (NLAM) . GT . 0945 . 0 ) WR I TE ( 6  *  40 34 ) 

IF (P (N LAM ) . LE. 0945 . 0 • AND.P ( NL AM) . GT . 0915. 0 ) W  R I T  E (6*4  0  0  5 ) 

IF (P( NLAM) .LE. 091 5. O.AND.P (NLAM) .GT .0885.0 ) WRITE ( 6  *  4036 > 

IF (P (NLAM) .LE. 0885.0. AND.P (NLAM) . GT  .OR  5 5 .0 ) WR I TE ( 6  *  403 7 ) 

IF (P (NLAM) .LE. 085 5. O.AND.P (NLAM) .GT .0825 .0 ) WR ITE ( 6  *  4038 ) 

IF (P (NLAM) .LE. 08  25 .O.AND.P ( NLAM) .GT .0795 .0 ) WR ITE ( 6  *4039 ) 

IF (P( NLAM) .LE. 0795 .O.AND.P (NLAM) . GT .0765 . 0 ) WR I TE ( 6  *4040) 

IF  ( PC  NLA  M)  .  L  E  •  0765  .  O.AND.P  ("NL  AM  )  .  GT.  073  5 . 0  )  WR  I  TE  (6*4041  ) 

IF (P (NLAM) .LE.0735 .O.AND.P (NLAM) .GT. 0705.0) WRITE (6*4042 ) 

IF (P (NLAM) .LE.0705 .O.AND.P (NLAM) .GT • 06  7  5 • 0 ) WR I TF ( 6  *  404  3 ) 

IF (P( NLAM ) .LE. 0675 .O.AND.P (NLAM) .GT .0645 .0 ) WR I TE ( 6  *  4044 ) 

IF (P( NLAM) .LE. 0645. O.AND.P (NLAM) .GT.  0615. 0) WRITE (6*4045 ) 

IF (P( NLAM) .LE.0615. O.AND.P (NLAM) .GT. 0685.0) WRITE (6 *4046 ) 

IF ( P (NLAM) .LE. 0585.0. AND.P ( NLAM ) .GT .0555 .0 ) WR ITE ( 6  *  4047 ) 

IF (P (NLAM) .LE. 0555.0. AND.P (NLAM) .GT.0525.0) WRITE (6 *4048 ) 

IF (P (NLAM) .LE. 0525 . O.AND.P (NLAM) .GT .0495 . 0 ) WR I TE ( 6 *4049 ) 

IF (P (NLAM) .LE. 0495 .O.AND.P (NLAM) .GT.0465 .0 ) WR ITE ( 6 *4050 ) 

IF (P (NLAM) .LE. 0466.0. AND.P (NLAM) .GT .0435 .0 ) WR ITE ( 6  *  405 1 ) 

IF (P (NLAM) .LE. 0435.0. AND.P (NLAM) .GT .0405 .0 ) WR I TE ( 6 *4052 ) 

T E ( P (NLA M) .LE. 0405 • 0. AND. P ( NLAM). G T . 03 7  5.0) WR IT F (6 *406 3 ) 

IF (P (NLAM) .LE. 0376.0. AND.P (NLAM) .GT. 0345.0) WRITE (6*4064 ) 

IF (P( NLAM) .LE. 0345 .O.AND.P ( NLAM) . GT . 03 1 5 . 0 ) WR I TE ( 6  *  405 5 ) 

IF (P( NLAM) .LE. 03 15. O.AND.P (NLAM) . GT . 0285 .0 ) WR I TE ( 6  *  4056 ) 

IF (P (NLAM) .LE.0285. O.AND.P (NLAM) .GT.0255.0) WRITE (6 *4057) 

IF (P( NLAM) .LE. 025  5.0. AND.P (NLAM) .GT .0225 .0 ) WR I TF ( 6  *  405 8 ) 

IF (P (NLAM) .L E • 02  2 5 .0. AND.P ( NLAM ) .GT. 0 1 9 5. 0 ) WRIT E (6  *40 59l 
IF (P (NLAM) .LE. 01 95.0. AND.P (NLAM) . GT . 0 1 6 5 . 0 ) WR I TE ( 6  *  406 0 ) 

IF (P (NLAM) .LE. 01 65.0. AND.P (NLAM) .GT .0 1 3  5 .0 ) WR T TF ( 6  *  406 1 > 

IF ( P (NLAM ) .LE. 01 36. O.AND.P (NLAM) . GT .01 05 .0 ) WR I TF ( 6 *406 2 ) 
IF(P(NLAM) .LE. 0105.0. AND.P (NLAM) .GT. 0075.0) WRITE (6 *4063 ) 

IF (P (NLAM) .LE. 0075 .O.AND.P ( NLAM) .GT .0045 .0 ) WR I TE ( 6  * 4064 ) 

IF  (P(NLAM)  .LE.004  6 .0.  AN  D  .  P  (  NLAM  )  .GT  .0  o  1  5 .0  )  WR  I  TF  "(6  *  A  0  6  5 ") 

IF (P (NLAM) .LE. 0015 .0. AND.P (NLAM) . GT • ( -001 5 • 0 > ) WR I TF ( 6  * 4066 > 

IF (P( NLAM) .LE. ( -00 1 6 . 0 ) . AND. P ( NL AM ) .GT. (-0046.0)  ) WR I TE ( 6  *  4067 ) 
IF (P (NLAM) .LE. ( -0045 . 0 ) .AND . P ( NLAM ) .GT . ( -0076 . 0 ) ) WR I TE ( 6 ♦ 4061 ) 
IF (P (NLAM) .LE. ( -0075 . 0 ) .AND . P ( NLAM ) .GT • ( -0105 • 0 )  ) WR I TE ( 6  * 4069 ) 
IF (P( NLAM) .LE. ( -0105 • 0 ) .AND.P ( NLAM ) .GT. (-0135.Q)  ) WR I TE ( 6  *  407 0 ) 
I F ( P ( NL AM ) .LE. ( -0135 . 0 ) .AND . P  ( NLAM ) .GT . ( -0166. OT) WRITE T 6  *  4071 T 
IF (P (NLAM) .LE. ( -01 65 . 0 ) • AND • P ( NL AM ) .GT • ( -0 1 95  #0 )  ) WR I TE ( 6  * 4072 ) 
IF (P (NLAM) .LE. ( -01 96 . 0 ) . AND. P ( NL AM ) .GT. (-02  26.0)  ) WR I TE ( 6  *  407  3  ) 
IF (P (NLAM) .LE. (-02 25.0) .AND.P (NLAM) .GT. ?-0255*0> ) WRI TE  <6  * 4074 ) 
IF (P (NLAM) .LE. ( -0255 . 0 ). AND. P ( NLAM ) .GT. (-0286.0) ) WR ITE ( 6  *  4075 ) 
IE (P (NLAM) .LE. (-0285 .0) .AND.P (NLAM) .GT. (-0315.0)  ) WR I TE (6  * 4076 ) 
TFTP'imTRH  . L E .  r-7)^TTT0TTCTTD7 M ( NL A'M'V VST - 0  3 4570 )  IWKTtE  (6*4077  ) 
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IF ( P ( NLAM ) .LE. ( -0345 . 0 ) . AND. P ( NL AM ) .GT. (-0375.0)  )  WR I TE ( 6  *  407  8 ) 
IF (P (NLAM) .LE. (“0878.0) .AND.P(NLAM) .GT. (-0408.0) ) WR I TE ( 6 ♦ 4079 ) 
IF (P (NLAM) .LE.  (-0408.0) . AND. P ( NLAM ) .GT. (-0488 >Q)  ) WR I TF ( 6  * 40 B 0 ) 
IF (P (NLAM) .LE. (-0488 .0) • AND • P (NLAM) .GT. (-0468 .0> ) WR ITE ( 6*4081  ) 

IF (P(NLAM) .LE. ( “0468 • 0 ) • AND.P ( NLAM ) .GT. (-0498.0) ) WR I TE ( 6  *  408 2 ) 
IF (P (NLAM) .LE. (-0495.0) .AND.P(NLAM) . GT . ( -08 2 8 . 0 )  ) WR I TE ( 6  *  40 8 8 ) 
IF (P (NLAM) .LE. (-0828.0) .AND.P(NLaM) .GT. (-0888.0) ) WR I TF ( 8 * 4u8 4 ) 
I F ( P ( NLAM ) .LE. ( -0888 • 0 ) .AND. P ( NLAM ) .GT. (-0888.0) ) WR I TE ( 6  *  408 8 ) 
IF (P( NLAM) .LE.  ( -0888 • 0 ) . AND . P ( NLAM ) .GT • ( -06 1 8 • 0 )  ) WR I TE ( 6  *  408 6 ) 
IF(P  '(NLAM')' .LE.  (-0618.0)  •  AND.  P  ( NLAM")  •  G  T  •  (-06  4  8 . 0  )  )  W  R  I T  E  (  6%  4087) 

IF (P (NLAM) .LE. (-0648.0) . AND. P ( NLAM ) .GT . (-0678.0) )WRITE (6*4088) 
IF (P (NLAM) .LE. (-0678.0) .AND.P (NLAM) .GT. (-0708.0)  ) WR I TE ( 6 * 40 R 9  ) 
IF (P (NLAM) .LE. (-0708.0) .AND.P(NLAM) .GT • (-0735*0 >  * WRITE (6 #4090) 
IF (P( NLAM) .LE. (-0788.0) • AND. P ( NLAM ) .GT • (-0768.0)  ) WRITE (6*4091  ) 
IF (P (NLAM) .LE. ( -0768 . 0 ) . AND. P ( NLAM ) .GT. (-0798.0)  ) WR I TE ( 6  *  409 2 ) 
I F ( P ( NLAM ) . LE . ( -0798.0 ) . AND • P (NLAM ) .GT . ( -0828 . 0 )  ) WR I TE ( 6  *  409  3 ) 
IF (P( NLAM) .LE. ( -0828 . 0 ) . AND. P ( NLAM ) .GT. (-0888.0) ) WR I TE ( 6  *  4094 ) 
IF (P (NLAM) .LE. (-0868.0) .AND.P(NLAM) .GT. (-0888.0) ) WR I TE ( 6 ♦ 409 8 ) 
IF (P (NLAM) .LF. ( -0888 . 0 ). AND.P ( NLAM ) .GT. <-0915 *6 > ) WR I TF ( 6  *  40 9 6 ) 
IF  (P  (NLAM)  .LE.  (-0918.0)  .AND.P  (NLA.M)  .GT.  (-0948.0)  )  WR  I  TE  (  6  *  4097  ) 
IF (P (NLAM) .LE. ( -0948 . 0 ) .AND . P ( NLAM ) .GT. (-0978.0) ) WR I TE ( 6  * 4098 ) 
IF ( P ( NLAM) .LE. (-0978.0 ) .AND.P ( NLAM ) .GT . ( -1005 .0 ) ) WRITE (6*4099 ) 
IF (P (NLAM) .LE. ( -1008 . 0 ). AND. P ( NLAM ) .GT . ( -1 03 8 . 0 )  ) WR I TE ( 6  * 4 1 00 ) 
IF (P(NLAM) .LE. (-1038.0) .AND.P(NLaM) .GT. (-1068.0) )WRITE (6*4101) 
IF  (P  (NLAM)  .LF.  ( -1068 . 0  )  .AND.P  (  NLA.M  )  .GT.  (-1098.0)  )  WR  I  TE  (  6  *  4 1 0  2  ) 
IF (P( NLAM) .LE. ( -1095 . 0 ) .AND.P ( NLAM ) .GT. ( -1 1 25 • 0 ) ) WR I TE ( 6 *4 1 03 ) 
IF (P( NLAM) .LE. ( -1 1 25 • 0 ) . AND. P ( NL AM ) .GT . ( -1 1 5 5 • 0 )  ) WR I TE ( 6  * 4 1 04 ) 
~ IF ( P ( N  LAM) . LE. ( - 1 1 55.0 ) . A  ND. P(NL A M ) . G T . ( - 1 185 .0> )WRI TF (6*4105 ) 

IF (P (NLAM) .LE. (-1185.0) .AND.P (NLAM) .GT. (-1215.0)  ) WR I TE ( 6  *  4 1 06 ) 
IF (P(NLAM) .LE. (-1 ?1 8.0) .AND.P (NLAM ) .GT. (-1245.0)  ) WR I TE ( 6  * 4 1 07 ) 
IF (P( NLAM) .LE. (-1245.0) .AND.P (NLAM) .GT. (-1278.0)  ) WR I TE ( 6  * 4 1 08 ) 
IF(P (NLAM) .LE. ( -1275 • 0 ) .AND. P ( NLAM ) .GT. (-1305.0) ) WR I TE ( 6  *  4 1 09 ) 
IF (P (NLAM) .LE. ( -1305 .0 ) .AND.P (NLAM ) .GT • ( -1335 • 0 ) ) WR I TE ( 6  *  41 1 0 ) 
IF (P(NLAM) .L E . ( - 133  5.0 ). AND.P ( NLAM ) . G  T. (-1365.0) >WR I TE ( 6*41 11 ) 

IF(P (NLAM) .LE. ( “1 365 • 0 ). AND • P ( NLAM ) .GT • (-1395.0) ) WR I TE ( 6  *  4 1 1 2 ) 
IF(P (NLAM) .LE. ( -1395 • 0 ). AND.P (NLAM ) .GT • ( -142 8 .0 > ) WR I TE ( 6  *  41 1 3 ) 
IF (P(NLAM) .LE. (-1425.0) .AND.P (NLAM) .GT. ( -1 455 .0 )  ) WR I TF ( 6  *  4 1 1 4 ) 
IF  (P  (NLAM)  .LE.  ( -1455 .0  )  .AND.P  (  NLA.M  )  .GT  •  ( -1485 .0  >  )  WR  I  TE  (  6  *  4 1 1 5  ) 
IF  (P(  NLAM)  .LE.  ( -1485 . 0  >.  AND  •  P  (NLA.M  )  .GT  •  (-1515.0))  WRITE  (6*4116) 
I FTP < NLAM) .LE • (-15 16 . 0 ) .AND.P ( N  LAM) . GT . ( - 1 545 .0  TTW RITE (6*4 1 17 ) 
IF (P (NLAM) .LE. (-1545.0) .AND.P (NLAM) .GT. ( -1575 .0 )  ) WR I TE ( 6  *  41 1 8 ) 
IF(P (NLAM) .LE. ( “16  76 . 0 ) . AND . P ( NLAM ) .GT . (-1605.0)  ) WR I TE ( 6  *  4 1 1 9 ) 
IF (P( NLAM) .LE. ( -1 606 • 0 ). AND. P ( NL  AM ) .GT. ( -1 63  6 . 0 )  > WR I TE ( 6  *  4 1 2 0 ) 
IF (P( NLAM) .LE. ( -16 35 . 0 ) . AND . P ( NL AM ) .GT . ( -1665 . 0 > ) WR I TE ( 6  *  4 1 2 1 ) 
IF (P (NLAM) .LE. (-1665.0 )• AND.P (NLAM) .GT . ( -1 695 .0 )  ) WR I TE ( 6  *  4 1 2 2  ) 
IF (P ( NLAM) .LF. ( -1695.0 ) .  AND . P ( NL  AM ) .GT .  (  -1725 .0 ) ) WR I TE ( 6  *  4 i 2 3 ) 

IF(P(NLAM) .LE. ( -1 726 . 0 ) .AND.P ( NLAM ) .GT. ( -1758 . 0 )  ) WR I TE ( 6  *  41 24 ) 
IF (P (NLAM) .LE. ( -1755 . 0 ) • AND . P ( NL AM ) .GT . ( -1 78  6 .0 )  )WR I TE ( 6 ♦ 4125  ) 
IF (P( NLAM) .LE. (-1785 .0) .AND.P (NLAM) .GT • ( -1 8 1 8 . 0 ) ) WR I TF ( 6  *  4 1 26 ) 
IF  (P  ('NLAM)  .GT.  1818  .O.OR.P(  NLAM)  .LT.  (-1815.0)  )  WR  I  TE  (  6  * 4 200  ) 

4006  FORMAT(005X*1H. ) 

4  007  FORM  AT ( 0  0  6  X  *  1 H  .  ) 

4008  FORMAT (007X*1H.) 

4009  FORMAT (008X*1H.) 

4010  FORMAT (009X*1H. ) 

4011  FORMAT (010X*1H.) 

4012  FORMAT  ( 01 IX  *  1H.  )  __  _  _  _ _ __ 
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4014  FORMAT t  0 13X 
4018  FORMAT (  0 1 4X 
4016  FORMAT (  0 1  5X 


4017  FORMAT ( 016X 

4018  FORMAT ( 0 1 7X 

4019  FORMAT ( 0 18X 

4020  FORMAT ( 0 19X 

4021  FORMAT ( 020X 

4022  FORMAT ( 021 X 
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1H.  ) 
1H.  ) 
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1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
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4024 

4025 
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FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 


022X 

023X 
0  24X 
025X 
0  26X 
027X 
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1H.  ) 
1H.  ) 
1 H  •  ) 
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1H.  ) 
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FORMAT ( 0  28X 

FORMAT ( 0  29X 
FORMAT ( 030X 
FORMAT ( 0  3 1 X 
FORMAT ( 032X 
FORMAT ( 033X 
FORMAT ( 0  34X 
FORMAT ( 035X 
FORMAT ( 036X 
FORMAT ( 037X 
FORMAT ( 038X 
FORMAT ( 039X 
FORMAT (040 X 
FORMAT ( 041 X 
FORMAT ( 042X 
FORMAT ( 043X 
FORMAT (  044X 
FORMAT ( 045X 
FORMATT046X 
FORMAT ( 047X 
FORMAT ( 048X 
FORMAT ( 049X 
FORMAT ( 050X 
FORMAT ( 051X 


1H.  ) 

1H  •  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1HTT 
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1H.  ) 
1 H  •  ) 
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1H. 

1H. 

1H. 

1H. 

1H. 
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FORMAT (0 52X 
FORMAT ( 0  53X 
FORMATt  054X 
FORMAT ( 055X 
FORMATt 056X 
FORMAT ( 057X 
FORMAT TO 5 8X 
FORMATt  059X 
FORMAT (060X 
FORMAT  t  061X 
FORMAT ( 062X 
FORMAT  t  063X 
FORMAT t064X 
FORMAT ( 065X 
FORMATt  066X 
FORMAT  t  067X 
FORMAT  f 068X 
FORMATt 069X 
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( •  #  X8I0 )TAMP 
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FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT  (_ 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 


07 1  X 
072X 
0  73X 
074X 
075X 
076X 
077X 
078X 
079X 
08  OX 
0  8 1 X 
0  82X 
083X 
0  84X 
085X 


1H.  ) 
1H.  ) 
1H.  ) 


1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1  H  •  ) 
1H.  ) 
1H.  ) 
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4094 
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4097 
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4104 
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FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
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FORMAT ( 
FORMAT ( 


086X 

087X 
088X 
089X 
090X 
091  X 
092X 
093X 
094X 
095X 
096X 
097X 
098X 
0  99X 
100X 
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10?X 
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) 
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) 
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1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1  H  •  ) 
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4105 

4106 

4107 

4108 

4109 
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FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 


1  04X 
1 05X 
106X 
107X 
1  08X 
1 09X 


1 1  OX 

111X 
11?X 
113X 
1 14X 
115X 


1H.  ) 

1H.  ) 
1H.  ) 
1H.) 
1H.  ) 
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4113 

4114 

4115 
41  16 


FORMATS 

FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 


4117 

4118 

4119 

4120 

4121 
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format ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 


116X 
1 17X 
118X 
119X 
120X 
1  21  X 


1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
1H.  ) 
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4124 

4125 

4126 
4200 


FORMAT ( 

FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
RETURN 
TNT5 - - 


1  2?X 

1  23X 
124X 


1  H  •  ) 

1H#  ) 
1H.  ) 


1 25X ♦ 1 H« ) 

1  OX  »  24H  VALUE 


TOO  LARGE  TO  PLOT) 


(  •  Hr*Xr’"r  )  T  '  C'  r  ,\ 

(  •  H  r  •  y  rTO  )  TAMOf)^  rrr.s 

'  ,Mr,xrrn)T^orn 
(  ,*■■  r  9  XA'“  r  )  T  A»*aqq  pr"  ,\ 
(.Hr,X^  a 

(,Hr.Y^^)T4'*goi  r^PA 
(  .  ‘J  r .  yrrr  )  TM/gn^  or-  *\ 
(.Hr*  XR^r  )  TAMqni  o^n  a 
(  .Hr.  XPVO)  TftMOT-1  OP  A  A 

( .H T . XPR° ) T AMOqP  rpo.N 

(  .h  r .  x  r^r  )  tavao^  ^rpa 
(•Hr.x^njr.Mam  pr^a 
( ,Hr»YfooiTA  »ao i  ar  p  a 

( . H r • XAQ p ) T Awqm  PP^A 
{  .Hr^X^riTAMPOl  *R^A 
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( .  H r# XVRP ) TAMqn^  PRO  A 
(•^r.XRRr)TAVQfn  P  p  n  A 
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SIBETC  ATLED 

SUBROUTINE  DELTA 
C 

C  THIS  SUBROUTINE  CALCULATES  THE  SHEAR  FORCES  BETWEEN  EACH  GAUGE 

C  LINE.  SUBROUTINE  FORCE  PRECEEDS  SUBROUTINE  DELTA.  I  AND  IN  MUST 

C  BE  D EEINED  DIREC TLY  R EEORE  C A LL  D E L T A . 

C 

C  REQUIRED  IN  COMMON  ST ( 30 *20  *  5 )  .DEL ( 24 ♦ ?0 ♦ 5 ) , DELT ( 24  * 20 * 5 ) 

C 

COMMON  DR(30*20*5) * ST (30*20*5) *DEL( 24.20*3 ) >DELT (24*20*5 > , 

1H0RC (30*5) *HORT ( ^0  *5 )  *SHOR ( 30.5 )  ,RhOR (30*5),TEIMO(30*5),EIMO(20)> 
1EXM0 (30*5) >SMOM(30*5>  »RMOM(30.5) *EC ( 20 >  *ET ( 20 ) *P C 20 ) *  A ( 20T»Y ( 20 ) , 
1X(30)*BB1(30,5)*RR2(30»5)*BR3(30,5)*C](30,6),C2(30*8>,C3(20,5)9 
1 AA1 *AA2 ♦ AA3  *AA4>AA5  * AP2*AP3  *AP4> AP5  *RP1 *BP2  >RP3  * 

1SH2 * SH3 *  SH4 , SH5 . SM2  9 SM3  *  $M4 , $M5  *  tH1 , JH2  * TH3 , TH4 , TH5 , TH6 , TH7 , TMT V 
1TM2*TM3*TM4*TM5*TM6*TM7*DIV»B1*B2*B3*DET  ♦ 

1B*D>CAL  *TL *GAGL*LIMI T .LOLIM  * L AML  * LSP AN » LOAD *PLOAD *  SPAN  * NGL *NLAM, 
1NL  *NNL  *  NGLM  *NGLG  *NGLL  *NGLR  *R ♦ I ,TN ,NNN ,NORUN“ 

C 

DO  Sno  NL  =  2  *NNL 
DO  502  NGL  = I *NGLL  *  2 
DO  503  NLAM=1*20 
IE (NORUN. LE. 2 ) GO  TO  507 
IF ( NORON. GE. 3) GO  TO  509“ 

507  DEL(NGL,NLAM,NL )=-(ST (NGL-2 *NL AM ,NL ) -ST ( NGL ,NLAM ,NL ) ) *1000.0 
DEL (NGL*NLAM*NL)=DEL(NGL*NLAMfNL) /(B*GAGL) 

DELT (NGL  »NL AM *  NL  )  =  0 • 0 
GO  TO  503 

5  09  DEL (NGL,NLAM,NL )=- (ST (NGL-4  >NLAM,NL )-ST ( NGL,NLAM,NL ) ) *1000.0 
DEL ( NGL ♦ NLAM  *  NL ) = OF  L T N  GLT N  LAM , N L )  77  B*2.0*GAGL ) 

DELT (NGL.NLAM,NL)=0.0 

503  CONTINUE 
502  CONTINUE 
500  CONTINUE 

1  =  1  +  1 

DO  504  NL=2VNNL 

DO  506  NGL= I  * NGLR  * 2 
DO  505  NLAM=1,20 
I E ( NOR UN .LE • 2 ) GO  TO  513 
IF (NORUN. GE. 3 )G0  TO  515 

513  DEL (NGL ,NLAM,NL ) =- ( ST (NGL-2  *NLAM ,NL )-ST ( NGL ,NLAM ,NL ) ) *1000.0 
DEL  ING LVNLA M  * NITI^DELTNG L7NL4M  ,N  L77TD  *G A G L ) 

DELT (NGL  *NL AM , NL ) =0 • 0 
GO  TO  505 

515  DEL (NGL  *NLAM,NL ) =- ( ST ( NGL-4  *NLAM ,NL  > -ST ( NGL ,NLAM ,NL )  ) *1000.0 
DEL (NGL  ,NLAM,NL ) =DEL ( NGL ,NLAM,NL ) / (B*2.0*GAGL ) 

DELT (NGL  *NLAM  »NL )=0.0 

505  CONTINUE- 

506  CONTINUE 

504  CONTINUE 
1=1-1 

DO  516  NL=2*NNL 
DO  518  NGL=T*NGLG 

- DO  5?U  NLAM=TT»TO - -  " 
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DEL (NGL .NLAM.NL ) =  -DEL (NGL.NLAM.NL ) 

520  CONTINUE 
518  CONTINUE 
516  CONTINUE 

DO  508  NL  =  2  *NNL 
DO  510  NGL= I *NGLG 
NLAM*1 

DELT(NGL*NLAM,NL)=DEL (NGL.NLAM.NL) 

DO  51?  NLAM=?,10 

D  F  L  T ( N  G  L  *  N  LA M  .  N  L  )  =  DEL  T ( N G L  *  N  L  A  M - 1 .  NT ) +DEL (NGL.NLAM.NL) 

51?  CONTINUE 
NL AM=?0 

DFLT(NGL.NLAM.NL)=DEL (NGL.NLAM.NL) 

514  CONTINUE 

NL AM=NL AM-1 

DELT ( NGL  * NLAM*  NL  ) =  DEL  T ( NGL *NL AM+l *NL  ) +DEL ( NGL .NLAM.NL ) 
IFfNLAM.GT.il )G0  TO  514 
510  CONTINUE 
508  CONTINUE 
RETURN 
END 


SIRFTC  CFRROC 

SUBROUTINE  CORREC 
C 

C  THIS  SUBROUTINE  ELIMINATES  THE  DISCONTINUITY  AT  THE  MID  DEPTH  OF 

C  THE  BEAM.  SUBROUTINE  DELTA  PRECEEDS  SUBROUTINE  CORREC.  TOTAL 

C  SHEARS  ARE  CORRECTED  SHEARS. 

C 

C  REQUIRFD  IN  COMMON  Y ( 20 )  * SHOR (30.5) *DELT (24*20*5) 

C 

COMMON  DR (30*20*5) » ST (30.20.5) *DEL ( 24  * 2 0  * 5  ) *DFLT ( 24  *  20  *  5 >  ♦ 

1H0RC (30*5) *HORT ( 30*5 ) *SHOR ( 30*5 ) ,RHOR( 30*5 )*TEIMO(30*5)*E!MO(?0)* 
1FXM0 (30*5 ) *  SMOM ( 30*5  >  *RMOM ( 30*5 )  *EC ( 20 ) *FT( 20 ) *  P (2  0 ) *  A (20)  *  Y ( 20 ) * 

IX (30) ,BB1 (30*5) *RB2 (30.5 ) *BB3 (30.5 ),C1(30*5)*C2(30,5)*C3(30*5), 

1AA1 .AA2  *AA3.AA4*AA5  *AP2*AP3*AP4*AP5  *BP1 *BP2  *RP3  »  _ 

1 SH2 . SH3  * SH4 . SH5 . SM2 . SM3 . SM4 . SM5 . TH1 ,TH2 , TH3 , TH4 , TH5 , TH6 , TH7 , TM 1 , 
1TM2.TM3 ,TM4*TM5  * TM6 * TM7 * D I  V  * B1 *B2*B3*DET, 

1B*D*CAL*TL  *GAGL  *LIMIT  *LOLIM*LAML,LSPANJ,LOAD.PLOAp  *  SPAN  *  NGL  *  NLAM . 

1 NL , NNL , NGLM ,NGLG ,NGLL »  NGLR , R , I , IN , NNN . NORuN 
C 

_ _ NLAM«1 

Y(NLAM) =D/?.0-TL/2.n 
DO  550  NLAM  =  2  *  9 
Y(NLAM) =Y(NLAM-1 )-TL 
550  CONTINUE 

NL AM=1 0 

Y ( NLAM )=Y(NL AM-1 )-0.75*TL 
NL  AM  =  1 1 

Y(NLAM) =TL/4.0 
NL AM= 1 2 

- V  ( NLAM  )  =V '( NLAM-1 )  4- ft  .7^'*Tr - - —  “ 


rf -o 


(!r'.  '  j )  i  J'--  r  .  va  j  - .  j  • )  i 

q  uj*  T  t  i^O ">  nC  ? 

!HMTTH03  Pfs 
> r  J 

JHH,  ct  jh  i  rn 
a  i n h •  I  —  idv  nr p  on 
r=M/\  jm 

f  M,  -A  j-.  if'M)  J  ■  1  =  (  l/,v/\  JH,  JDM)  T  H 

nr#c=^AJH  c  r  p  00 

M JH* f-MA JH* Je>M)TJ30»( JH#Ma JH* J3HITJ 

hjhithoo  cr? 

AC5ma JM 

(  JH.V'A  !H,J0H)  JlP,r(  JH,”A  JH,  J0H)T  J10 

HlHlTWO')  A  I  ^ 
r-MAJH»MA JM 

f  !  ■*#’  '  M  .  Jr )  J  -  r+  (  JM,  r+V*  v  *  J  0  M  )  T  J 1 1  =  (  JM*”A  JH  •  !  n  H  )  T  Jin 

AT*  0T  Of)  (  r  r.TD.MAJMJIT 

It  IMI  TMOO  nf? 
I'MITMOO  ?np 

(1IO 


oopqnn 

M5I9M  qMITUOqqiP 


lTqqr  JM  1HT  TA  YT  1 1 IM  I  TMOl?  I  0  ^HT  ?  ^  T  A  H  I  M  I  J  1  ?M  I  TUOflflUr*  ’IH^ 
TOT  •  '  r  JO  )?  ATJ30  IMH  ■  #M  IT 

•  ppn  HP  fnr->qqaoA  ?1A  ?9A1H? 

(  ’P^AC'T  nnrP.nnqoM.Mnov  mom  vim  HT 

.  (  jqn,  (  a.rc,Of  )T?,  (^,A^pf  )fln 

'  r  >  y/  T  T#  (  P,  O  r  )  C'*T  IT,  f  p.  or  )  80‘Jfl.  (  a .  or  )  qo»-p  .  (a,Of)  T  qOH,  f  ?  ,  n  r  )  nqOM  r 

A  r  )  V.  (  >  ,  )  a  .  (  n<; )  q,  (  nc)  Tl .  (  nr  )  01,  f  q.  nr  )  MCH'P#  (  ’,nr]  MOM? ,  (  ?  .  nr  )  omyi  r 

.  (  ’  #  n  r  )  c  0  .  (  >,AH  O,  (  )  H,  I  n. ''f  )  tQ^*  (  p»  n^)  cqo,  (  P,  off  )  TAP.  pf  JXT 

•  r  i,  r  q3f  f  q n  •  P  q  A  *  A 1 A  *  P  q  A  •  q  *W  ?  A  A  •  A  A  A  ,  C  A  A  *  c;  A  A  »  r  A  A  J 

•  .  IT  *  AHT • £HT • SHT • f HT  *  3M? • AM? ♦ #  $M?  *  £H  .  ■ 

r.P.VTO.T^T,  }VT  .  3*'T,  AVT#  PVT,  ^’t; 

iM*J£)M*MA  >*  \0QA  .  ■  It  JMAJtMIJOJtTTMl  J»  JOAO#  J  t  '  •  • 

M»r  ■)  -I .  m  v*.  1  j »  t  ,  .  '  r  I  J-iH.  -1  JO'/.  -*•*)/.,  JMM,  '/r 


rr*«A  JH 


3 

*> 

0 

0 

3 

0 


3 


n*S\JT«  .  '  :  |M  J 

c#c  =  ^ajh  n a nn 
J  T  —  (  r-MA  JH  J  Ys  r  WA  JM  )  V 

IjiHItmo-)  np  a 


n  r  =  M  A  j >/ 


JT^T.n-f  f  -MA  JM  )  Y-  (»AJW)Y 


f  fsMA JH 
O.A\ JT= (YA  JM)  V 
9  r-tWA  JH 

|T*ar  .  r>  +  (  r-‘M  i /My -('//•  ]!/■) 


C-32 


DO  652  NLAM= 1 3*20 
Y(NLAM) =Y ( NL AM- 1 ) +TL 

552  CONTINUE 

DO  5  54  NGL  =  ] ,NG L M 

DO  556  NL=1 ,NNL 
5HOR ( NGL ,NL )=0.0 
556  CONTINUE 
554  CONTINUE 

DO  553  NLAM=1,?0 

Y ( NLAM )  =  ( (D/2.0 ) -Y ( NL  AM ) ) / (D/2.0 ) 

553  CONTINUE 

DO  558  NGL= I  * NGLG 
DO  560  NL=2 ,NNL 
NLAM=1 0 

IF(DELT(NGL,NLAM,NL).GT.DELT(NGL,NLAM+l*NL) )GO  TO  566 
IF  ( DELT (  NGL  »NL AM »NL ) • EQ.DEL T ( NGL .NLAM+ 1  * NL ) ) GO  TO  564 

IF (DELT (NGL,NLAM,NL) • LT.DELT (NGL.NLAM+1 ,  NL ) )GO  TO  562 
566  SHOR(NGL,NL )=ABS( DELT (NGL, NLAM. NL) -DELT (NGL, NLAM+1 ,NL) ) /2.0 
DO  568  NLAM=1,10 

DFLT(NGL,NLAM,NL)=DELT(NGL,NLAM,NL)-(SHOR(NGL,NL)*Y(NLAM) ) 
568  CONTINUE 

DO  570  NLAM=1 1,20 

DFLT (NGL, NLAM, NL )=DELT (NGL, NLAM, NL )+(SHOR( NGL, NL )*Y( NLAM) ) 
570  CONTINUE 
GO  TO  564 

56?  SHOR (NGL ,NL  )  =  ABS ( DELT ( NGL ,NLAM,NL)-DELT( NGL ,NLAM+1 , NL ) ) /2.0 
DO  572  NL AM= 1 ,10 

DELT(NGL,NLAM,NL“)=DELTTNGL,NLAM,NL  )  +  (SHOR(NGL  ,NL  )  *  Y  ( NL  AM )  ) 
572  CONTINUE 

DO  574  NLAM=1 1,20 

DEL T ( NGL, NL AM, NL )= DELT (NGL, NL AM, NL)-( SHOR ( NGL, NL )*Y( NLAM) ) 
574  CONTINUE 
564  CONTINUE 
560  CONTINUE 

558  CONTINUE 
RETURN 
END 
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